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Abstract 

The quest for the Higgs boson of the Standard Model, which was a cornerstone in the physics 
programme at particle colliders operating at the energy frontier for several decades, is the subject 
of this review. After reviewing the formulation of electroweak symmetry breaking via the Higgs 
mechanism within the Standard Model, the phenomenology of the Higgs boson at colliders and the 
theoretical and phenomenological constraints on the Standard Model Higgs sector are discussed. 
General remarks on experimental searches and the methodology of statistical interpretation are 
followed by a description of the phenomenology of Higgs-boson production and the corresponding 
precise predictions. The strategies of the experimental searches and their findings are discussed for 
the Large Electron Positron Collider (LEP) at CERN, the proton-antiproton collider Tevatron at 
Fermilab, and the proton-proton Large Hadron Collider (LHC) at CERN. The article concludes 
with the description of the observation of a Higgsdike boson at the LHC. 
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1 Introduction 



Electroweak history in a nutshell 

More than 40 years ago the standard theory of electroweak interaction p] , the Glashow-Salam- Weinberg 
model, emerged from a longer theoretical development, starting from the old Fermi theory in the 
1930s, being generalized to a phenomenological model of intermediate massive vector bosons, and 
finally receiving the form of a spontaneously broken gauge theory in the 1960s. In those days the last 
bottleneck in the mathematically consistent construction of the electroweak theory was the merge of the 
successful Yang-Mills structure for the interaction of vector bosons with fermions with the experimental 
fact that the weak vector bosons possess mass. A historical milestone in this development was laid with 
the realization that a spontaneous breakdown of the gauge symmetry, driven by a self-interacting scalar 
field, can lend mass to gauge bosons — an idea nowadays known as the "Englert-Brout-Higgs-Guralnik- 
Hagen-Kibble (EBHGHK) mechanism" or simply "Higgs mechanism". 

In the early 1960's it was realized by Nambu, Goldstone, Salam and Weinberg [2] that the sponta- 
neous breakdown of a global continuous symmetry necessarily leads to the existence of a massless scalar 
("Goldstone") particle — a statement known as "Goldstone's theorem". This fact prevented theorists 
from associating the postulated massive weak vector bosons with a broken symmetry, because no cor- 
responding Goldstone bosons had been observed. Later it was realized by Brout, Englert j3], Higgs [I], 
and Guralnik, Hagen, Kibble [5] that a spontaneous breakdown of a local continuous symmetry does 
not require Goldstone bosons, but rather their degrees of freedom deliver the longitudinal polarization 
modes of the gauge bosons that become massive. The understanding of this phenomenon, called the 
EBHGHK or Higgs mechanism, lead to the mathematical formulation of the class of spontaneously 
broken gauge field theories. In his second paper [I] PW. Higgs stated "It is worth noting that an essen- 
tial feature of the type of theory which has been described in this note is the prediction of incomplete 
multiplets of scalar and vector bosons." It seems that it is this statement on the particle structure of 
such models that finally gave the Higgs particle its name. 

The qualitative picture of the Higgs mechanism is that all massive elementary particles^ interact 
with the vacuum, which is characterized by a non-vanishing expectation value of the Higgs field. This 
interaction with the vacuum prevents massive particles from acquiring the speed of light, in analogy 
to light rays proceeding through matter. At the same time, the Higgs field itself admits quantum- 
mechanical excitations which appear as new scalar particle types, the Higgs bosons. 

Based on the understanding of the Higgs mechanism to lend masses to gauge bosons, Weinberg and 
Salam successfully completed Glashow's model for the unified electroweak interaction to the sponta- 
neously broken SU(2)xU(l) gauge theory that still represents the heart of the electroweak part of our 
Standard Model (SM) of particle physics today. The remaining part of the SM is the gauge theory 
of strong interactions, known as quantum chromodynamics [6j, which is based on the unbroken SU(3) 
colour gauge symmetry. The final breakthrough of gauge theories in particle physics came in the early 
1970's with the proof of their renormalizability, both for unbroken and broken symmetries by 't Hooft, 
Veltman [7J, Lee and Zinn- Justin [8|. This step raised gauge theories generally to the level of mathe- 
matically consistent quantum field theories, a fact that serves as our basis to work out predictions for 
collider experiments at a level of precision that is essentially only limited by our technical capabilities to 
evaluate higher orders in perturbation theory. Of course, some non-perturbative input is needed, such 
as parton distribution functions for hadronic collisions, but conceptually there is a solid field-theoretical 
basis. 



1 Here we understand that "elementary" means "point-like" and thus non-composite particles, because the mass of the 
latter typically receives large contributions from their binding energy. In the case of nucleons, the binding energy even 
comprises the major part of the mass. 
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Prom theory to Higgs phenomenology and searches 

The particle content of the SM comprises three charged leptons, three corresponding neutrinos, and six 
quarks in the sector of fermions. These "matter fermions" interact via the exchange of gauge bosons — 
the photon 7 for electromagnetic interaction, the weak gauge bosons Z and W , and the eight gluons 
of the strong interaction. Apart from this particle content, which is experimentally established after 
the discoveries of the top quark at the Fermilab Tevatron in 1995 and of the r neutrino at the Fermilab 
DONUT experiment in 2000, there is just one postulated electrically neutral Higgs boson in the SM 
which is the phenomenological footprint of the Higgs mechanism. Since masses play the role of coupling 
strengths between particles and the vacuum, and Higgs bosons appear as "vacuum excitations", the 
coupling of the Higgs boson to any other particle is predicted to be proportional to the particle's mass 
term in the field equations. A phenomenological confirmation of the Higgs mechanism, thus, requires 
finding the Higgs boson(s), measuring its/their quantum numbers such as spin and electrical charge, 
and finally verifying the proportionality of the Higgs coupling to the mass of the attached particle. 
Once a value for the hypothetical Higgs-boson mass is assumed, the profile of the SM Higgs boson 
is completely fixed, and precision predictions are possible for production and decay rates and signal 
process kinematics. Based on these precise predictions experimental searches have been performed for 
the signature of the SM Higgs boson at colliders operating at the energy frontier. It took 25 years 
after the formulation of the Higgs mechanism until a significant mass range could be probed with the 
start of the operation of the Large Electron Positron Collider (LEP) at CERN in 1989. The search was 
continued at the Tevatron proton-antiproton collider from 2002 to 2011 at Fermilab. In 2010 the Large 
Hadron Collider (LHC) started to take data in proton-proton collisions at unprecedented CM energies 
with a primary goal of finally answering the question of whether a SM Higgs boson is realized in nature 
or not. 



2 The Standard Model Higgs boson 

Before we enter a thorough discussion of the phenomenology of the Higgs boson in the SM of particle 
physics, we briefly review the general structure of the model. Here we focus on the features of the 
SM that are most relevant for Higgs-boson phenomenology and refer to standard textbooks such as 
Refs. [9TTT3] for more details and theoretical background. 



2.1 The gauge structure of the Standard Model 

The Lagrangian Csm of the SM can be divided into four different parts, 

£SM = £YM + £ferm + + ^Yuk, (1) 

which will be discussed one by one in the following. We start with the first two, which contain the 
dynamics of the gauge bosons and matter fermions as predicted by the unbroken gauge symmetry, 
before we turn to the latter two which involve the Higgs field. 

The Yang-Mills part £ym describes the genuine dynamics of the gauge fields, i.e. their free propa- 
gation as well as their self-interactions, 

£ YM = -~W; V W^ - \b, u B^ - -Gl v G^\ (2) 



where 





= w 






= d,B u 


— d v B^, 




= d,G a v 


— d v G a ^ 



h Ji k — 1) 2, 3, 



fahc G\Gl, a,b,c =!,...,&, (3) 
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are the field-strength tensors of the gauge fields for the SU(2)i group of the weak isospin 1^, £? M for 
the U(l)y of weak hypercharge Y w , and G a ^ for the SU(3) C of colour. The respective gauge couplings 
of these groups are denoted g, g', and g s , and the structure constants of the non-abelian groups SU(2) 
and SU(3) are e l ^ k and f abc , following the usual notation. 

The interaction of the gauge fields with the fermions is encoded in 

£ ferm = M L p* L + $e R fl^t R + MqP*q + iVv#>Vw + % R , (4) 

where L = (ve L ,£h) T are the left-handed SU(2)i doublets of charged leptons I — e, r and neutrinos 
vi = u e , z/ M , u T , Q = (ml, c?l) T the left-handed SU(2)i doublets of up-type quarks u — u, c, t and down- 
type quarks d = d, s, b, and £r, %, g?r are the respective right-handed SU(2)i singletso The actual 
interaction is contained in the covariant derivative 

= d, + igP w W; + ig'Y w B, + ig^G^ (5) 

where J^, Y w , and T c a are the generators of the respective gauge groups in the representation of the 
fermions they act on, i.e. = cr*/2 (a 1 = Pauli matrices) for the left-handed SU(2)j doublets and 
14 = for the right-handed singlets, the weak hypercharge Y w (= numbers) is related to the relative 
electric charge Q by the Gell-Mann-Nishijima relation Q = I^ + Y w /2, and T c a = A a /2 (A a = Gell-Mann 
matrices) for SU(3) C quark triplets and T c a = for the leptons. The requirement that the coupling 
structure of the photon is parity blind and proportional to QipAip, as in quantum electrodynamics, 
identifies the photon field and the Z-boson field by the rotation 



-Afj, J V s w c w J \ B ^ 



(6) 



with the weak mixing angle 9 W and electric unit charge e fixed by 



9 99' 



cos6> w = c w = yjl - 4 = =, e = . (7) 

lg 2 + (g') 2 Jg 2 + (g') 2 



Once the photon is identified among the gauge fields, it is easy to see that the fields 

W^ = (W^iW 2 )/V2 (8) 

correspond to the charged weak gauge bosons of charge ±e. 

Note that the part of the SM described so far, encoded in £ Y m + £ferm, does not involve any mass 
terms. Naive gauge-boson mass terms, such as W^W 1 ^ for the W bosons, obviously violate gauge 
invariance. Since left- and right-handed fermions transform differently under SU(2)ixU(l)y gauge 
transformations, naive fermion mass terms oc (ipj h tpf R + V'/rV'/l) f° r a fermion / are also ruled out by 
gauge invariance. The introduction of particle masses, while still maintaining the gauge invariance of 
the dynamics of the model, requires extensions of the theory. 



2.2 Electroweak symmetry breaking and the Higgs mechanism 

The Higgs part 

£ H = (D^{D»<S>) - V($) (9) 

of the Lagrangian extends the experimentally well established particle content of the SM by the complex 
scalar SU(2)i doublet $ = (</> + , </>°) T of weak hypercharge Y w ^ = 1, so that <f) + carries charge +e and 

2 Right-handed neutrinos are omitted, since they do not play any role at high-energy colliders. 
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</>° is neutral. In total $ involves four real degrees of freedom. The self-interaction of $ is described by 
the potential 

V($) = -At 2 ^) + ^($ f $) 2 , (10) 

whose form is constrained by gauge invariance and renormalizability of the model. The latter is guaran- 
teed by the polynomial structure of degree four in the field components, the former by the $ dependence 
via the combination leaving only the two real free parameters fi 2 and A in V. While A > is re- 
quired by vacuum stability, the sign of /x 2 is deliberately taken positive in order to force a non-vanishing 
vacuum expectation value (vev) $o of Minimizing V yields the condition 

*J*o = y, v = 2^£. (11) 

Requiring that the vev is electrically neutral, forces the upper component of $o to vanish, i.e. <3>o is fixed 
up to a phase, with the usual choice <3>o = (0, vj V2) T . This freedom in choosing the vev $o of $ reflects 
the spontaneous breakdown of the SU(2)iXU(l)y symmetry down to the remaining electromagnetic 
U(l) em invariance. Splitting off the vev from <3>, 



®-{<f> = (v + H + ix)/V2)' (12) 

we reparametrize $ in terms of the real physical Higgs field H and the unphysical would-be Goldstone 
boson fields <fr + and x> which are complex and real, respectively. The fact that <p + and x do not 
correspond to physical states can already be seen by the fact that they are connected to the vev by 
gauge transformations; in fact one can always find a gauge, known as the "unitary gauge" , in which + 
and x vanish. Making use of this gauge and inserting the parametrization ( fl2|) of $ and the covariant 
derivative (jSJ) with = <t % /2, Y Wj <f> = 1, T c a $ = into the Higgs Lagrangian (J9j), we find 



which in particular contains bilinear terms in the gauge fields W^, Z and in the Higgs field H, i.e. mass 
terms for the corresponding weak gauge bosons W ± and Z as well as for the Higgs boson H. Identifying 
these masses according to 

M w = ^, M z = ^, M H = v/V, (14) 

we can eliminate the parameters /i 2 , A, v completely and get 

Ai,u- g au g c = \{dHy-\MlH 2 + M^W+W->»+ l -MlZ»Z» 

+ gMwHW+W-^ + —H 2 W+W-'" + 9 —HZ„Z^ + ^-H 2 Z.Z^ 

4 2c w 4c4 

4M W 32M^ ' K J 

where we have not spelled out an irrelevant constant. In summary, the SM makes the following important 
phenomenological predictions in the Higgs sector which can be tested experimentally: 

• Associated to the Higgs field H, a physical neutral, spinless particle of mass Mh is postulated — the 
Higgs boson. Since it can be viewed as some kind of vacuum excitation, it carries the quantum 
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numbers of the vacuum and is thus even with respect to CP symmetry. Note that Mr is the only 
free SM parameter that is tied to a property of the Higgs boson, while the other parameters are 
fixed by the weak-gauge-boson masses and the gauge couplings. Theoretical and phenomenological 
constraints of Mr are discussed below. 

• The ratio p = M^j / (c^M|) |14^l is equal to one, which is a non-trivial relation among the weak- 
gauge-boson masses and the gauge couplings g = e/s w and g' = e/c w . While p = 1 in this form 
is used to define c w = Mw/Mz in the process of renormalization to all orders, it nevertheless 
has important phenomenological consequences. For instance, on top of the Z-boson resonance, as 
measured at LEP and SLD via e + e~ — > Z — > ff, an effective weak mixing angle can be defined 
for each fermion species /, usually quantified via sin e f eS , which can be measured from various 
asymmetries. The ratios pf = MyV/ (cos 2 0{$M%) thus are predicted to be equal to one up to 
radiative corrections, a fact that is experimentally confirmed at a level of better than 10~ 3 [T5] . 

The lowest-order property p = 1 is not shared by all possible scalar sectors that can be employed 
to lend masses to the W ± and Z bosons. If not accidental or forced by fine-tuned parameters, it 
is rather a consequence of a "custodial symmetry" [T6] of the scalar sector. In the SM this is an 
SO (4) symmetry of £h with respect to the exchange of the four real components of $, which holds 
up to U(l)y gauge interactions and differences in the fermion masses within SU(2)i doublets, both 
inducing custodial-symmetry-breaking effects, however, only in higher orders. 

• The model predicts couplings of the Higgs boson to a massive weak gauge boson V = W, Z pro- 
portional to My/v, because the HV'V couplings originate from the factor (1 + H/v) 2 multiplying 
the gauge-boson mass term in £ H - Owing to the square in this factor, quartic couplings of two 
Higgs bosons and two gauge bosons are also predicted, which are proportional to My/v 2 . 

• Finally, triple and quartic Higgs-boson self- interactions are predicted, both scaling with My. Since 
these couplings are in one-to-one correspondence with the shape of the Higgs potential that drives 
electroweak (EW) symmetry breaking, an experimental reconstruction of these couplings from an 
analysis of scattering processes would be part of an ultimate phenomenological confirmation of 
the Higgs mechanism. However, since the relevant processes involve multi-Higgs-boson final states 
with very low cross sections, the LHC will at best be able to give qualitative results here. 

2.3 Yukawa couplings and fermion masses 

Renormalizability and gauge invariance of the SM Lagrangian allow for so-called Yukawa couplings of 
the Higgs doublet $ to all fermions. The most general form of these interactions is 

Cyuk = -^ L GWv R $ - ^ Q G u ip u ^ - * Q G^dR$ + h.c, (16) 

where "h.c." means hermitian conjugate and $ = i<7 2 $* = ((0°)*, — 4>~) T denotes the charge-conjugate 
Higgs doublet with quantum numbers opposite to $. The matrices Gf (/ = £,u,d) represent arbitrary 
complex 3x3 matrices, i.e. at first sight Cy n k involves a large number of free parameters. However, most 
of them turn out to be not physically relevant and can be transformed to canonical values or eliminated 
by appropriate field redefinitions. We first note that each term in £yuk involves terms that are bilinear 
in the fermion fields because $ and $ contain a constant piece in the form of the vev v. More precisely, 
the non-diagonal elements of Gf mix the left- and right-handed parts of the different generations of 
fermion type / — £,u, d, where £ generically stands for charged leptons, u for up-type quarks, and d 
for down- type quarks. Owing to these mixing terms, a fermion of flavour /j (i = 1,2,3) of the ith 
generation oscillates into fj of the other generations (j ^ i) even during a free propagation in space and 

3 The ratio p is called /3 in Ref. [H]. 



7 



time. This oscillation can be removed upon transforming the existing "flavour basis" (iftf Tl ,ipf T2 , ^Us) 
of left- and right-handed fields (r = L, R) into a "mass basis" (^1,^2,^3) with a unitary matrix 

U, 

j>fr,i = U{fipf T j, f = *,u,d, r = L,R, (17) 
where in this process the matrices Gf receive a diagonal form, 

U*-G f {Uf*)i = ^diag{m fl ,m h ,m f3 }. (18) 

The diagonal value my\, which can be chosen non-negative by convention, is the mass of the fermion 
/j. For completeness we mention that the left-handed neutrino fields are transformed with the same 
unitary matrix as their charged counterparts; this is possible as long as we work in the approximation 
of three massless (i.e. mass-degenerate) neutrino's. The effect of this field redefinition on the whole 
SM Lagrangian can be summarized easily: (i) The coupling matrices Gf are replaced by their diagonal 
form (fl8|) ; (ii) all fermion fields if>f Ti are replaced by their counterparts ipf ri of the mass basis; (iii) 
the only remnant of the U matrices is the appearance of the Cabibbo-Kobayashi-Maskawa matrix 
V = U UL (U dL y in fermion chains of type ^ nL . . . V^d L and of in . . . V^ UL . In other words, only 
charged-current interactions receive modifications by V, while neutral currents remain unchanged. In 
the following we adopt the common convention to omit the clumsy hats on fermionic fields and assume 
the use of the mass basis. 

In the unitary gauge, where the would-be Goldstone fields are absent, the Yukawa Lagrangian takes 
the simple final form 

1 + — j, (19) 

where the sum over / runs over all fermion flavours of all generations. This form shows a distinctive 
footprint of the Higgs mechanism in the fermionic sector: The Higgs boson couples to each fermion 
/ of mass rrif with the strength yf = irtf/v. Moreover, the coupling is the one of a pure scalar, i.e. 
the coupling to fermions does not have any pseudo-scalar admixture proportional to 75. Testing these 
features offers a possibility to empirically tell a potential Higgs candidate from scalar particles predicted 
by other models. Alternative models with non-minimal Higgs sectors often predict new pseudoscalars as 
well, or even scalar particles without definite CP quantum numbers. Moreover, the strict proportionality 
of the Yukawa coupling strength to the fermion masses might be broken, as it is for instance the case 
in (type-II) Higgs doublet models, where the proportionality factor between yj and my is different for 
up- and down-type fermions. 



2.4 From input parameters to predictions for collider experiments 
Running parameters and renormalization scale 

Like any gauge-field theory, the SM involves ultraviolet (UV) divergences that are removed in the 
process of renormalization which ties the input parameters of the theory to measurable quantities 
at some renormalization scale /zr. The input parameters, which are taken from experiment, in this 
sense depend on the arbitrary value of /zr as well as on more details that fix the actual procedure, 
thereby defining the renormalization scheme. The fact that predicted physical observables, on the 
other hand, cannot depend on /ir is expressed in terms of renormalization group equations (RGE) for 
each scheme (if it allows for a flexible renormalization scale). Solving these RGE for any field-theoretical 
quantity naturally involves the concept of running, i.e. /ir- dependent, parameters, such as the famous 
strong coupling q; s (/xr). The running of each parameter from one scale to another can be predicted 
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in perturbation theory. In actual predictions of observables, thus, there is an explicit and an implicit 
dependence on //r, which compensate each other order by order in perturbation theory — the former 
resulting from the renormalization of the considered observable at the scale /zr, the latter from the 
input values matched to renormalized quantities at /ir. A residual /xr dependence arises from orders of 
perturbation theory that are not completely taken into account. If an observable involves a typical scale 
Q, such as a scattering cross section at energy Q, it is advisable to fix //r in the vicinity of Q, since this 
minimizes potentially large missing corrections that involve powers of logarithms ln(Q//iR). If a process 
involves many different scales, the benefit of adjusting /iR is limited. For the precise mathematical 
formulation of this concept we refer to standard literature on quantum field theory such as Refs. [9HT3]. 

The input parameters of the SM 

As is obvious from the construction of the model, the free input parameters of the SM are the gauge 
couplings g, g', g s , the parameters /x 2 and A of the Higgs sector, the fermion masses nif, and the CKM 
matrix V. For phenomenology it is much more convenient to take instead the following parameters as 
input: the electromagnetic coupling a = e 2 /(47r), the strong couplings constant a s = g 2 /(47r), the weak 
gauge-boson masses M w and M z , the Higgs-boson mass M H , and finally rrif and V. The masses can all 
be defined as pole masses, defined from the locations of the particle poles in the respective propagators, 
but for the heavy quarks it is often useful to switch to a running mass at some appropriate scale. For 
Higgs physics, the CKM matrix plays a minor role. The couplings a and a s , however, have to be chosen 
thoughtfully. The strong coupling is usually defined as running coupling « s (/xr) in the so-called MS 
scheme. Very often the value of a s at the Z pole, a s (Mz), is used as numerical input and transferred to 
some other value /ir using the RGE. 

For the electromagnetic coupling a basically the choice is between three different values: the fine- 
structure constant a(0) ~ 1/137, the effective value a(M z ) « 1/129, where a(0) is evolved via RGE 
from zero-momentum transfer to the Z pole, and an effective value derived from the Fermi constant 
leading to a Gfl = a/2G^M^(1 - M^/Mf)/^ « 1/132, defining the so-called "G M -scheme" . The various 
values of a differ by 2—6%. Very often the actual value of a can be adjusted in such a way that large 
universal EW corrections are already absorbed into the lowest-order prediction. Even different values 
of a can be appropriate in one calculation, but care has to be taken that the same choice is taken 
within gauge-invariant subsets of diagrams. An external photon (virtuality Q 2 = 0) always effectively 
couples with a(0), while internal photons with virtuality Q 2 at a high energy scale effectively couple 
with az(Q), which is much closer to a(Mz). On the other hand, the couplings of W and Z bosons 
should be parametrized with , which does not only take into account the running from Q = to the 
EW scale, but also universal effects from the p-parameter. Following these rules, in particular, avoids 
perturbative instabilities due to the appearance of light-quark masses, as e.g. discussed in Ref. [IT] . 

Predictions for hadronic collisions and factorization scale 

Hadronic collisions are treated within the QCD-improved partem model, whose full description is beyond 
the scope of this review (see e.g. the standard textbooks [TlTri3] ). In this framework hadronic cross 
sections are calculated in the factorized form 



&hihi — 




where the parton distribution functions (PDFs) fai/hj^h Vf) are generalized probability densities for 
finding a parton d{ (quarks, antiquarks, gluons, photons) with momentum pi = XiPi in the hadron hi 
which has momentum P«. The PDFs comprise the soft physics, which is dominated by small energy scales 
and not accessible by perturbation theory, and are process independent. The partonic cross sections 
da aia2 , on the other hand, contain the full information on the hard scattering reaction of the partons 
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into a specific final state and can be calculated within perturbation theory. The separation of the soft 
and hard domains is performed at the factorization scale /Xp, which effectively is the upper bound on 
transverse momenta up to which outgoing hadronic particles are considered as parts of the struck hadron. 
Similar to /xr, the scale /xf is arbitrary to a large extent, and predictions — if calculated to arbitrary 
precision — cannot depend on it. In d<r aia2 (xiPi, X2P2, A*f) the value of /xp can be trivially changed, in 
the PDFs the change of /xf is governed by the so-called Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) evolution, which ensures the /XF-independence of observables. Again a truncation of the 
perturbative series for da aia2 , however, leads to some residual /xp dependence of observables, which is 
part of the intrinsic uncertainty in predictions. 

The PDFs cannot be calculated from first principles, but are fitted to a large variety of data re- 
sulting from lepton-hadron and hadron-hadron collisions. State-of-the-art PDF sets are provided by 
several groups: CT10 [18J, MSTW2008 [19], NNPDF2.1 [20J , ABKM11 ]21j, GJR08 [22], and HERA- 
PDF1.5 [23J. Generally, there is good agreement between predictions, e.g. for the LHC, based on those 
PDF sets, in particular between CT10, MSTW2008, and NNPDF2.1, and differences between the PDF 
sets are continuously analyzed. More details and references on PDFs and their application to Higgs 
physics at the LHC can be found in Refs. [23H2B] . In practice, each perturbative order requires its own 
set of PDFs for consistency, i.e. a LO cross section <tlo has to be evaluated with LO PDFs and the 
corresponding one-loop running a s (/x), onlo with NLO PDFs and two-loop a s (/x), etc. In this context, 
QCD corrections are often quantified by the K-factor 



resulting from the cross sections consistently calculated in different orders. While up-to-date PDF sets 
are available in LO, NLO, and NNLO QCD, no current PDF set includes EW corrections, i.e. EW 
corrections are both missing in the PDF fit to data as well as in the /xp-evolution of the PDFs. The 
older PDF set MRST2004QED [27] includes these effects, but is outdated otherwise. It is, thus, more 
advisable to use up-to-date PDFs, since ignoring EW PDF corrections typically introduces uncertainties 



Parametric and theoretical uncertainties 

Parametric uncertainties arise from the errors of the input parameters of the SM model. For observables 
that involve strongly interacting objects the largest errors usually originate from the uncertainties of the 
strong coupling constant a s and of the masses of the heavy quarks that are possibly directly involved. 
In observables based on the EW interaction often the uncertainty of the electromagnetic coupling a 
becomes significant if it is derived from the value of the running coupling a(Q 2 ) at high energies, such 
as a(M|). Parametric errors of different observables are usually correlated, which has to be respected 
when their effect is combined. 

For hadronic collisions the PDFs are very often the largest source of parametric uncertainties, which 
are also correlated with a s . The various groups who deliver regular updates of PDFs fitted to data 
provide different "error sets" of PDFs whose spread quantifies the PDF uncertainty. These error PDFs 
also allow for a combined assessment of the PDF+a s error of observables. A proposal how to combine 
errors from PDFs of different groups can be found in the report of the PDF4LHC Working Group [24]. 

Theoretical uncertainties can be — to some extent — quantified by the "scale uncertainties" mentioned 
above. The compensation between explicit and implicit /xr and/or xxp dependences remains necessarily 
incomplete in perturbative calculations, because in practice not all (but only few) perturbative orders 
can be calculated completely. Thus, a residual scale dependence remains in predictions, and, if the scales 

4 Though outdated, the MRST2004QED PDF set is, however, the only one that provides a photon PDF and, thus, 
offers the opportunity to at least approximately calculate the effect from partonic channels involving initial-state photons. 



K 



0"(N)NLO 



(21) 



0"LO 




10 



are varied within a reasonable range, the corresponding variation in the result quantifies theoretical 
uncertainties originating from missing higher-order effects. One should, however, keep in mind that 
not all missing corrections are sensitive to the residual scale dependences. For instance, missing finite 
corrections depend on only very weakly, but might be large. Thus, overdoing a minimization of the 
residual scale dependence does not remove the uncertainties, but rather tends to hide them. 

Finally, parametric uncertainties and theoretical uncertainties have to be combined. Since theoretical 
uncertainties have no true statistical meaning that is safely tied to a probability distribution, the two 
errors should not be added in quadrature, but linearly. A possible procedure to do that for Higgs 
observables at the LHC was, e.g., proposed in Ref. [25] . 



2.5 Theoretical and phenomenological constraints on the Higgs sector 

Although the mass Mh of the Higgs boson represents a free parameter of the SM, it cannot assume 
arbitrary values. In the following we briefly summarize the bounds on Mh from above and below that 
originate from theoretical considerations and phenomenological constraints. 



The heavy-Higgs limit of the SM 

Before looking more closely into the phenomenology of Higgs bosons, it is interesting to discuss the 
behaviour of the model in the limit Mh — » oo, which means that Mh is supposed to be much larger 
than any other mass or energy scale involved in the considered physical system. In this limit the Higgs 
boson is effectively removed from the physical particle spectrum, but its vev still delivers the masses to 
the weak gauge bosons and massive fermions. Formally the Higgs potential tends to a delta function 
proportional to — $0^0); i- e - the absolute value of the Higgs doublet field is constrained to 

its vev. This constrained scalar sector is more conveniently described in the framework of non-linear 
representations of the doublet $, where the physical Higgs boson plays the role of the radial excitation 
in the space of its four real components and gets fixed by the delta function. 

To illustrate this, we first rephrase the Higgs Lagrangian OH]) in matrix notation, 

£ H = hr{(D^)\D^)} - V($), = + i|aW;$ - i^ 3 5 M , (22) 



where the hermitian matrix field $ is given by 



• = -!. 



v 



+ H-i X )/V2 cj> + 



7# + F)1 + 10V K -0 (v + H + i X )/V2) = ^)> (23) 

with a 1 again denoting the Pauli matrices and the real scalar fields <p l , which are defined by (f) ± = 
((p 2 ± if/) 1 )/ a/2 and x — —<fi 3 - In this linear representation the Lagrangian and the meaning of the scalar 
field components are exactly the same as described before, only the bookkeeping is somewhat different. 
The non-linear representation of <fr parametrizes the matrix field according to 

* = ±=(v + B)U, [/ = expj^j, (24) 

which results from the linear parametrization upon a non-trivial field transformation. In this parame- 
trization the field H is gauge invariant. Inserting the non-polynomial field <fr into the Higgs Lagrangian 
( 12~2"|) . leads to a theory that is physically equivalent to the SM with the more common linear Higgs 
representation [29], i.e. the S'-matrix is the same in the linear and non-linear representations, although 
the Green functions do not coincide. The new form of the Higgs Lagrangian is 

= \iv + H)hr{{D,U)\D»U)} + ±(d»H)(d»H) + ^(v + Hf - A(„ + H f. (25) 
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Note that here the fields 0* of the unphysical Goldstone bosons, which are contained in the matrix U, 
appear only in their dynamical part oc tr\DU\ 2 , but not in the Higgs potential, where only (y + H) 
survives owing to the unitarity of U. For this reason this form is very attractive for the discussion of the 
heavy- Higgs-boson limit, even in gauges different from the unitary gauge, which corresponds to U = 1. 

To lowest order of perturbation theory, the heavy- Higgs-boson limit of £h,iJ is reached upon setting 
the Higgs field H to zero, since all Higgs propagators are suppressed by factors 1/M| i n amplitudes — a 
suppression that cannot be compensated by the M§ factors in the Higgs self-couplings, since there are 
always more Higgs propagators than Higgs self-couplings in amplitudes without external Higgs bosons. 
In other words, a heavy Higgs boson decouples in tree-level amplitudes. The resulting theory without 
a Higgs field is called the gauged non-linear sigma model [30] and described by the scalar Lagrangian 

Cm = - A vhi{{D,U)\D^U)} = CuAH = 0). (26) 

This model is equivalent to the SM in the unitary gauge without Higgs boson, but still represents 
an SU(2)ixU(1)y gauge theory. The model, though, is physically unsatisfactory, since it is neither 
renormalizable, nor does its S-matrix respect unitarity. In fact these two short-comings are closely 
related, as shown in Ref. [3Tj . 

Beyond leading order, a heavy Higgs boson does not decouple anymore in amplitudes, a fact that 
is visible in the different UV behaviour of the SM with and without Higgs boson — the latter being the 
gauged non-linear sigma model. While the SM with Higgs boson is renormalizable, i.e. closed in the 
UV sector, the gauged non-linear sigma model is non- renormalizable, i.e. incomplete in the UV sector 
and thus at best can serve as an effective field theory up to some finite energy. Thus, removing the 
Higgs boson from the SM particle content upon increasing its mass more and more, Mh effectively 
works as UV cutoff scale for the UV divergences of the gauged non-linear sigma model. Specifically, 
using regularization in D dimensions the uncancelled one-loop UV poles 1/(D — 4) of the sigma model 
are in one-to-one correspondence with one-loop corrections involving In Mh in the heavy-Higgs SM. At 
one loop, these non-decoupling effects of a heavy Higgs boson have been quantified in terms of effective 
Lagrangians in the literature [32], keeping also finite terms that are not suppressed by powers of Mh. At 
higher loop orders, non-decoupling effects can even grow with powers of M H , as e.g. seen in heavy-Higgs 
corrections oc M H and M^ to the p-parameter at the two-loop [33] and three- loop level [31], respectively, 
suggesting a scaling law of M^ L ~ 1 ' at L-loop order, possibly modified by powers of lnM H . 

The SM is of course not the only possible "UV closure" of the non-renormalizable non-linear o~- 
model. Any renormalizable model containing it bears some dimensionful quantity A, such as a particle 
mass, that acts effectively as cut-off for UV divergences, like the Higgs-boson mass in the SM. In the 
limit of large A, the same structure of In A terms will, thus, show up as in the heavy-Higgs SM with 
Mh playing the role of A. 

Constraints on Mh from unitarity arguments 

The unitarity of the S-matrix implies bounds on the high-energy behaviour of partial waves of 2 — > 2 
scattering processes, but a naive power counting of energy factors in amplitudes with longitudinally 
polarized weak gauge bosons Vl = Wl,Zl suggests that these bounds do not hold. The leading 
high-energy term in the longitudinal polarization vector Ey L (k) of an external Vl with momentum 
which behaves ~ k^/My, introduces a dangerous energy factor which seems to spoil the power 
counting. Unitarity is restored by cancellations that are governed by gauge invariance in the form of 
Slavnov-Taylor identities, which relate the fc M contraction of an external V M leg of an amplitude to 
another amplitude in which Vl is replaced by the respective Goldstone-boson field 0v Schematically 
an amplitude A^v L ... behaves like 

My L ..=eyUk)^ k0 ~ Mv ^I$ L .„ =AW.., (27) 
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where Ai^... is the "scattering amplitude" of <py (defined with an appropriate phase) at high energies 
with a decent scaling law in the scattering energy, as demanded by unitarity. Equation (|27|) expresses 
the so-called Goldstone-boson equivalence theorem [HUES]; whose formulation to all orders requires 
some care |36j . 

It is interesting to note that the equivalence theorem can also be formulated within the gauged non- 
linear sigma model [57]. Since the non-polynomial structure of the Goldstone-boson Lagrangian ( 12 6 p 
involves arbitrarily high powers of derivatives (momentum factors in amplitudes), however, the power- 
counting of energy factors in (1271) does not lead to a decent high-energy behaviour of the amplitude 
.Mv L ...- Amplitudes Aiy L ... in fact potentially violate unitarity, both in the gauged non-linear sigma 
model and in the SM with a heavy Higgs boson, Mh — > oo. In the unitary gauge, where no Goldstone- 
boson fields are present, the same conclusion is drawn from the "bad high-energy behaviour" of the W 
and Z propagators. 

The violation of unitarity can already be observed in the SM with finite Mh if the Mr value is large 
enough. The larger Mh, the more delayed is the necessary unitarity cancellation for large k°. If Mh is 
too large, the amplitude starts to "explode" before the damping due to Higgs-boson exchange sets in. 
In 2 — > 2 particle scattering, the most delicate unitarity cancellations occur in longitudinal vector-boson 
scattering Vl Vl — > Vl Vl whose tree- level amplitudes respect unitary only if M H < 1 TeV (38] . Within 
pure perturbation theory it is not possible to give sharp unitarity limits on Mh, because the range 
around Mh ~ 1 TeV represents the transition region where perturbation theory runs out of control in 
the Higgs sector due to the huge self-couplings. 

Triviality and vacuum-stability bounds 

It is interesting to apply the concept of a running coupling to the Higgs-boson self-coupling A. The 
running of A(/ir) is mainly driven by loop corrections to the if "^interaction involving the Higgs-boson 
self-couplings or top quarks with the largest Yukawa coupling. Present-day evaluations of A(/zr) involve 
(at least) two-loop effects and RGE. The behaviour of the resulting running A(/zr) at large values of 
/ir contains information on the scale A up to which the SM can exist as a self-consistent quantum field 
theory: 

• For a large Higgs-boson mass the Higgs self-interactions drive A(/zr) to larger values with increasing 
/ir, reflecting the lack of asymptotic freedom. Above some scale the coupling enters its non- 
perturbative domain, a fact that is most obvious at one- loop order where A(/xr) even reaches 
a so-called Landau pole and diverges. At two-loop order the divergence is mitigated to a UV 
fixed-point, but at three- loop order a Landau pole is expected to appear again. More details 
on the running and matching issues can be found in Ref. [3S]. Of course, the transition to the 
non-perturbative domain is a continuous process, so that only a generic scale A for the onset of 
this region can be given. The non-perturbative domain and the Landau pole could be avoided if 
the exchange of new particles (not part of the SM) in the loop corrections to the if 4 interaction 
counterbalances the too rapid running induced by the graphs containing the Higgs self-couplings. 
Thus, A can be interpreted as a lower bound on the new-physics scale if the SM is supposed 
to stay in its perturbative regime. Being perturbative or non-perturbative is certainly not a 
physical criterion on a field theory. However, the appearance of a Landau pole does not only 
signal non-perturbativity, but also a mathematical inconsistency. Beyond the Landau pole the 
self-consistency of the field theory seems to require a scalar self-coupling A that vanishes at all, 
leading to a trivial theory. The bound on A to avoid the Landau pole is, thus, also called the 
triviality bound [10] • We are not aware of a field-theoretical proof of the triviality of scalar field 
theories with 4 interactions, such as the SM, but the triviality hypothesis is supported by non- 
perturbative results from lattice field theory (see Ref. [H] arid earlier references therein). 
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Figure 1: Perturbativity (triviality) upper bound and vacuum-stability lower bound on Mr as function 
of the scale A at which new physics has to appear (at the latest) in order to avoid the bounds (plot taken 
from Ref. [B]). The two different versions of the perturbativity bound reflects the intrinsic uncertainty 
of the limit, and the bands of the (meta) stability bounds reflect parametric uncertainties. 

• For a small Higgs-boson mass the running of A(//r) is dominated by the top-quark loops which 
drive A(/zr) even to negative values for large /ir. This fact has drastic consequences on the Higgs 
potential that determines the non-vanishing vev of the Higgs field. In lowest order the potential 
( JlOjl still shows the asymptotics oc Aif 4 , but A(/xr) < for large /ir indicates that higher-order 
effects bend the Higgs potential V over to negative values for large Higgs fields H, i.e. the potential 
is not bounded from below anymore. Thus, a stable vacuum does no longer exist, and the whole 
model becomes inconsistent. A thorough theoretical analysis of the vacuum structure goes beyond 
a pure running-coupling analysis and takes into account the leading static (for Higgs fields at zero 
momentum) all-order-resummed corrections to the Higgs potential V, leading to the so-called 
effective Higgs potential V e s |42j . For small Mr the SM can escape the vacuum- stability bound [13] 
only if new-physics effects at some scale A modify the vacuum structure determined by V e s in 
such a way that a stable vacuum exists. This results in a lower bound on Mr for a given value 
of A. The bound can be relaxed to some extent by assuming a metastable EW vacuum with a 
lifetime longer than the age of the Universe, where the vacuum decay may be mediated either by 
zero-temperature quantum fluctuations or thermal fluctuations. 

Figure [IJ which shows the results of the recent two-loop RGE analysis of Ref. jH] , illustrates the upper 
and lower bounds on the SM Higgs-boson mass resulting from the hypothesis that the SM Higgs sector 
should be perturbative (and thus escapes the triviality bound) with a stable or at least metastable 
vacuum up to a potential new-physics scale A. A SM Higgs boson with M H < 200 GeV admits the SM 
to escape the non-perturbativity/triviality bound up to scales A < 10 12 GeV, far above any energies 
reachable in collider experiments. On the other hand, the LEP exclusion limit Mr > 114.4 GeV pushes 
the vacuum (meta) stability limit for the onset of new physics already to A ~ 10 6 GeV, which is out 
of reach for colliders as well. For a Higgs boson within the mass window 130 GeV < Mr < 170 GeV 
the SM would be surely consistent up to scales as high as the Planck scale A Planck = 2 x 10 18 GeV, 
where gravity effects should become as large as the strong and EW interactions. Finally, interpreting 
the recent LHC discovery as a SM Higgs boson with Mr = 126 GeV still allows the SM to be valid up 
to the Planck scale if our universe resides in a metastable phase. 
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Figure 2: The "blueband plot" (left, taken from Ref. [15]) showing the result of a global x 2 fit of the 
SM to precision data, projected onto the Mh axis and comparison of constraints on the top-quark and 
W-boson masses (right, taken from Ref. [16]) resulting from direct measurements and a global SM fit. 



Constraints from precision data 

Experiments at LEP1, LEP2, SLC, and the Tevatron have performed a large variety of high-precision 
measurements, such as the determination of the W, Z, and top-quark masses and widths, measurements 
of cross sections and various asymmetries in the process e + e~ — >■ 7*/Z — > f f at the Z pole, etc. These 
measurements characterize the past two decades as the era of EW precision physics. Parametrizing SM 
predictions as described in Section EU thus, renders a fit of the complete SM to these data possible, 
which highly constrains the SM input parameters (with the exception of the CKM matrix, where 
dedicated observables from flavour physics are needed). Figure [2] shows, on the l.h.s., the status of March 
2012 of such a fit projected on the Higgs-boson mass, together with the exclusion limits from the direct 
searches to be discussed later in detail. The blue band indicates theoretical uncertainties from missing 
higher-order corrections in the (numerous) underlying precision calculations, which were condensed into 
the state-of-the-art codes Zfitter [17] and TopazO [48] (see also Ref. [M] for details and references). 
The best fit value together with its 68% CL limits for the Higgs-boson mass is Mh = 94+24 GeV, being 
compatible with the open mass window 122 GeV < M H < 127 GeV (status July 2012, see Section [83!) 
which is obtained from the 95% C.L. exclusion limits of the LHC from below and above, respectively. 
The r.h.s. of Figure [2] illustrates the SM overall fit, as obtained by the Gfitter [46^1 collaboration, by 
its projection into the m t -Mw plane and compares the preferred fit region with the directly measured 
values of the top-quark and W-boson masses. The shown SM prediction for Mw as function of mt 
and Mh (diagonal lines for some fixed Mh values) is obtained from the measured muon lifetime (often 
translated into the Fermi constant G^). All these constraints on m t and Mw are perfectly compatible 
with the recent observation of a Higgs-boson candidate of a mass around M H = 126 GeV, to be discussed 
in detail below. 



5 Similar results have been presented in Ref. [50] , 
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2.6 Higgs-boson decays 

The Higgs boson of the SM is unstable and predominantly decays into the heaviest particle-antiparticle 
pair that is kinematically possible depending on the available energy Mh. The search for the Higgs 
boson, which has to be reconstructed from its decay products in detectors, thus depends on the Higgs- 
boson mass, since the relevant final states strongly vary with Mh. Since the branching ratio BR(H — > 
X) = Th^x/Thi which is the probability for the Higgs boson to decay into the final state X, weights the 
relevant Higgs-boson production cross section for the exclusive channel H — > X, theory has to provide 
precise predictions for all significant Higgs-boson decay widths — irrespective of whether the channel 
could be detected or not — since they all enter the total decay width T H , which is the sum of all partial 
widths r H -*.x- To get the global picture we show the SM Higgs-boson branching ratios and the total 
decay width in Figures [3] and HI where we separate the 1ow-Mh range that was accessible at the time of 
the e + e~ collider LEP (Figure [3]) from the Mh range that is challenged by the hadron colliders Tevatron 
and LHC (Figure HI). The relevant leading-order Feynman diagrams for the various decay channels are 
depicted in Figure [51 As anticipated above, the highest branching ratios typically correspond to decays 
into the heaviest particle-antiparticle pair that are kinematically allowed for a given value of Mh. 
However, as far as predictions are concerned, the regions where new channels become significant are 
theoretically delicate. For heavy-quark pairs the threshold regions are particularly sensitive to the 
heavy-quark masses whose proper perturbative treatment within QCD is highly non-trivial. For the 
Higgs-boson decays into the weak-gauge-boson pairs WW and ZZ, even the regions near and below the 
WW/ZZ thresholds become important which implies that the decays of the W/Z bosons have to be 
included in the predictions. 

Owing to the proportionality of the Higgs couplings to the masses of the decay products, the total 
Higgs boson width Th stretches over many orders of magnitude, viz. from 10 -4 GeV for Mh < 10 GeV 
to 10 3 GeV for Mh ~ 1 TeV. Because of the finite energy resolution of detectors, the width of a light 
Higgs boson is too small to be resolved in invariant-mass distributions of its decay products. On the 
theoretical side this fact justifies the approach of dealing with stable, on-shell Higgs bosons in production 
processes and factorizing the subsequent decay from the production. For large invariant masses of the 
Higgs-boson decay products resulting from H — > VV (V = W, Z), i.e. for M vv > 2M V > M H , off-shell 
effects can be quite large, 0(10%), but this region can be excluded by dedicated cuts [53]. On the other 
hand, a heavy Higgs boson would show up as an extremely broad resonance, whose width becomes of 
the same size as its mass for Mh IS 1 TeV. Theoretically this situation is very challenging, because 
production and decay processes do not factorize anymore from each other. Instead, a proper treatment 
of the broad resonance has to deal with the signal, consisting of Higgs-boson production, propagation, 
and decay, background, comprising non-resonant diagrams with the same final state as the Higgs-boson 
decay, and of interference effects between signal and background at the same time. Even the proper 
field-theoretical definition of mass and width of a heavy Higgs boson to parametrize the resonance 
becomes subtle. For more details about these issues, which are still under investigation, we have to 
refer to the literature (see e.g. Refs. [251[2ni[MlES] and earlier references therein). In the following we 
put the emphasis on the Higgs-boson mass range 100 GeV < Mh < 200 GeV, which is favoured by the 
overall fit of the SM to precision data and the results of the direct searches. 

To deliver precise predictions for the Higgs-boson decay widths and branching ratios, a huge effort 
was made by many theorists. The decay channels that are most important for Tevatron and the LHC 
are: 

• H -> // (mainly / = r, b, t) [SSrEB] 

In describing a Higgs boson decaying into bottom (or even lighter) quarks it is essential to base the 
Yukawa coupling on the running quark mass at the relevant scale, which is set by the Higgs-boson 
mass. For instance, for Mh ^ 100 GeV the transition from the pole to the running bottom mass 
mb(Mn) reduces the H — > bb partial decay width by ~ 60% or more. Starting from this improved 
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Figure 3: Branching ratios of the SM Higgs boson (left) and total decay width (right) for Higgs-boson 
masses accessible at LEP and before, calculated with the program Hdecay 
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Figure 4: Branching ratios of the SM Higgs boson (left, taken from Refs. [26l[52]), with the band 
widths illustrating the parametric and theoretical uncertainties, and total decay width (right, taken 
from Ref. [25]) m the Higgs-boson mass range accessible by the LHC. 




Figure 5: Leading-order diagrams for the various SM Higgs-boson decay channels, where Q denotes any 
heavy quark. 
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LO prediction, the perturbative QCD series, which is known in NLO [58] and beyond that even 
up to NNNNLO [Sn], shows nice convergence with a small residual scale uncertainty of ~ 0.1%. 
Recently, the NNLO QCD corrections to H — > bb became available for fully differential observ- 
ables [61] as well. Generally, a proper treatment of the qq threshold region deserves particular 
attention [59] . 

For the decay into top quarks, the full mass dependence of the tt final state has to be included, and 
the issue of a running mass is not as pronounced as for the lighter quarks. The QCD corrections, 
which are available at NLO [58j[62] and NNLO [63], turn out to be moderate, but still ~ 5% at 
NNLO. The decay of the top quarks should be taken into account in predictions as well, not only 
for precision, but also since it would offer the possibility to analyze the CP properties of a possible 
heavy-Higgs-boson candidate (see Refs. [0HEH] anc ^ references therein). 

The theoretical structure of the NLO EW corrections to H — > ff is very similar for all 

fermions, and their size is typical of the order of a few percent. Beyond NLO, some QCD and EW 
corrections that are enhanced by factors of G^m^ are known [67J as well as the leading two-loop 
heavy-Higgs-boson effects oc G^M^ [68] that reflect the failure of perturbation theory in the TeV 
range for Mh, where they are as large as their NLO counterparts. 

• H -> gg [SHE] 

The NLO QCD (two-loop) corrections to the H — > gg decay were first calculated as an asymptotic 
expansion in Mn/m t [71], which is applicable only for Higgs-boson masses below the tt threshold 
(Mh < 2m t ), and later generalized to the full mass dependence [72] . These corrections are very 
large, ranging from 70% to 40% for M H = 100 GeV to 1 TeV. QCD effects beyond NLO [73], 
which are known in the form of an expansion in small values of Mh valid up to Mh < m t , add 
roughly another 20% to the partial decay width of a light Higgs boson to gluons. Electroweak 
NLO corrections are known with the full mass dependence and amount to ~ 5% [741175] . Mixed 
QCD-EW beyond NLO corrections were calculated in the heavy-top limit, i.e. for light Higgs 
bosons, but found to be very small [76] . 

• H ->■ 77/7Z [M1E71I7JI711ES1IZM3] 

Similar to the gluonic Higgs-boson decay, the NLO QCD corrections to the H — > 77 decay were first 
calculated in the large-top-quark-mass limit [78] and later with the full mass dependence [721179] . 
For small Mh these corrections amount to only < 5%, but for masses above the tt threshold 
they can be as large as 50—100% around Mh ~ 600 GeV, where the LO prediction receives some 
suppression due to a destructive interference between quark and W-boson loops. QCD effects 
beyond NLO are also known in the form of an expansion in small values of Mh [82]. NLO EW 
corrections, which were presented in Refs. [74, 75, 81J, turn out to be of the order of a few percent. 
In the actual calculations, particular care has to be taken in the vicinity of the tt and WW 
thresholds, where the finite decay widths of the top-quarks and W-bosons in the loop have to be 
taken into account [74 1 175 1 150"]. 

Because of its lower phenomenological importance, the predictions to the decay H — > 7Z are much 
less advanced than in the 77 case. The NLO QCD corrections, which turn out to be below 1%, 
can be found in Ref . |83j . Many leading higher-order corrections are very similar or even identical 
in the 77/7Z cases. 

• h -»■ ww/zz 4/ [651IB1HSZ] 

Since the Higgs-boson decay channels into weak-gauge-boson pairs are also phenomenologically 
relevant in the WW/ZZ threshold region and below, i.e. for M H < 2M V (V = W, Z), the W/Z 
decays have to be taken into account in precision calculations. Moreover, it was extensively 
discussed in the literature [84] that the kinematical details (energy, angular, and invariant-mass 
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Table 1: Estimated relative uncertainties of Higgs-boson branching ratios due to theoretical plus para- 
metric errors for low and intermediate Higgs-boson masses, as given by the LHC Higgs Cross Section 
Working Group [26j|52]. 



distributions) of the four-fermion final state offer the possibility to extract the spin and CP 
properties of a Higgs-boson candidate. Radiative corrections have first been calculated for on-shell 
W/Z bosons [65l|85], the results of which are only applicable for Mh > (2My + some GeV), and 
later including the full off-shell effects of the W/Z bosons [SB] in the framework of the complex- 
mass scheme [SSJ for treating the W/Z resonances. For H — > WW/ZZ the heavy-Higgs-boson 
limit is of particular interest, because the dominating decay into longitudinal W/Z bosons leads 
to the above-mentioned rise oc M H of the total decay width Th- The leading heavy-Higgs-boson 
corrections to r H -s.ww/zz are proportional to G^M^ and Gj;M H at the one- and two- loop level [87]. 
The fact that this two-loop term becomes as large as the one-loop correction for Mh ~ 1 TeV 
signals again the onset of the non-perturbative regime for a heavy Higgs boson. The leading two- 
loop effect, which is negligible for Mh < 400 GeV and amounts to ~ 4% for Mh ~ 700 GeV, can 
serve as an estimate for the theoretical uncertainty for large Mh. For 100 GeV < Mh < 200 GeV 
the radiative corrections to the partial widths range from some percent to < 12%, depending on 
the four-fermion final state and can be even much higher in distributions. 

More details can be found in reviews or reports [25 | [26 | |69 | [B9]. 

State-of-the-art predictions can be obtained upon combining results from the two public programmes 
Hdecay [51] and Prophecy4f [SB] : 

• Hdecay is designed to calculate all relevant Higgs-boson decay widths and branching ratios 
in the SM and its minimal supersymmetric extension. On the QCD side, all relevant known 
corrections are included, while the EW corrections, which are smaller for most of the channels, 
are approximated. Specifically, the EW (two-loop) corrections to H — )• gg/77 are reconstructed 
via a numerical grid obtained from the full results of Ref. [74"|I75]. and the corrections to WW/ZZ 
(considered as two-body decays) are fitted to the corrections provided by Prophecy4f. 

• Prophecy4f is a Monte Carlo event generator for all four-particle decays H — > WW/ZZ — > Af 
that takes into account all interferences and all off-shell effects of the intermediate W/Z bosons, 
so that the obtained results are valid above, near, and below the gauge-boson pair thresholds. As 
an event generator, Prophecy4f can simulate the four-body decays with its full kinematics with 
NLO QCD and EW corrections; for leptonic final states unweighted events can be produced as 
well. For large Mh the leading two-loop EW corrections from Ref. [57] are included as well. 

All results shown in Figure [BJ and all but the decay widths for H — > WW, ZZ in Figure H] are based 
on HDECAY, while the results on H — > WW, ZZ are obtained upon summing over the relevant four- 
fermion final states in H — > WW/ZZ — > 4/ whose decay widths are calculated with Prophecy4f. 
More technical details about the calculations for Figure H] can be found in Refs. J25J|2BJ[52] from which 
the plots are taken. 



Mh [GeV] H -> bb t+t 



cc gg 77 WW ZZ 



120 3% 6% 

150 4% 3% 

200 5% 3% 



12% 10% 5% 5% 5% 
10% 8% 2% 1% 1% 
10% 8% 2% <0.1% <0.1% 
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The bands in the prediction of the branching ratios shown in Figure H] reflect the combined un- 
certainties due to missing higher-order effects in the calculations and parametric uncertainties in the 
input parameters, as assessed by the LHC Higgs Cross Section Working Group • Table [T] briefly 

summarizes the estimated total uncertainties for Mh values 120 GeV, 150 GeV, and 200 GeV, which 
are most interesting for SM Higgs physics at Tevatron and the LHC. It is important to note that the 
uncertainties of the branching ratios result from an interplay of the uncertainties of the corresponding 
partial decay width and the normalization to the total width Th- The decay widths into hadronic 
final states (bb, cc, gg) possess relatively large uncertainties via the parametric uncertainties from the 
heavy-quark masses and the strong coupling. As long as Th is dominated by r H ^bE this uncertainty is 
mitigated in BR(H — > bb), but feeds into the branching ratios of the other channels, as, e.g., clearly 
visible for H WW/ZZ at M H = 120 GeV. The opposite is observed at M H = 200 GeV, where T H 
is dominated by r H ^ww/zz, which are very precisely known, i.e. in this case the uncertainties of the 
other branching ratios entirely reflect the one of the respective partial widths. A similar estimate of 
the uncertainties of the branching ratios has been presented in Ref. [20] in the Higgs-boson mass range 
of Mh = 100—200 GeV, where somewhat larger uncertainties are found, mainly due to larger assumed 
errors on the b- and c-quark masses. 

For very low and very high Higgs-boson masses the theoretical uncertainties of the branching ratios 
become larger and larger and hard to quantify. Figure [3] shows the branching ratios without error bands, 
mainly because the predictions in the vicinity of the qq thresholds, i.e. for Mh ^ 15 GeV, become 
more of qualitative nature with uncertainties of 100%. This theoretical shortcoming was balanced 
experimentally by inclusive search strategies in the 1ow-Mh range, where the individual decay channels 
of the Higgs boson did not play a role. For large Higgs-boson masses the uncertainty of predictions 
increases more and more when the non-perturbative regime of Mh ~ 1 TeV is approached. 

3 General considerations for Higgs-boson searches 
3.1 Challenges 

The challenge in searches for the Higgs boson is to develop a selection that allows to increase the a 
priori very small signal-to-background ratio (see Figure |6]) to a level of the order 1/1 to 1/10. In order 
to reach this goal, sophisticated selection strategies have been developed over the last decades at the 
experiments at LEP, Tevatron, and the LHC. In the beginning only simple rectangular cuts were used 
to separate the signal process from background processes. By now many search strategies make use of 
multivariate techniques such as likelihood ratios, artificial neural networks, or boosted decision trees. 

The combination of production mechanisms and decay modes of the Higgs boson and the particles 
that are eventually produced in association with the Higgs particle yields a plenitude of final-state 
topologies which can be used to search for Higgs-boson production. The following considerations have 
to be taken into account when choosing the most promising search channels: sufficient a priori signal rate 
determined by cross section and branching ratios, means to trigger on signal events with high efficiency, 
means to suppress the background to an acceptable level, means to estimate the uncertainties in the 
background prediction, and in the case of exclusion also the uncertainties in the signal prediction. 
It is desirable to be able to construct a final discriminating variable to distinguish the two alternative 
hypotheses of background-only and signal+background in addition to simple event counting. An obvious 
choice is the reconstructed mass of the Higgs-boson candidate events. Depending on the final-state 
particles this is not always possible. Alternatively or in addition other powerful observables can be 
used, e.g. the response of a multivariate event classifier. 

For the discovery of a new resonance, which means that the probability to observe the actual number 
of events is less than 2.87 x 10~ 7 under the background-only hypothesis, only the background prediction 
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Figure 6: Signal and background cross sections at leading order (produced with Pythia [9T1I92] ) at 
e + e~ colliders (left) and at hadron colliders (proton-antiproton and proton-proton colliders). 



and its uncertainty need to be known. Wrong assumptions on the shape of discriminating observables 
for the signal process will lead to a reduced sensitivity. Hypothetically one observed event can yield the 
claim of discovery given the search is basically background free (~ 3 x 10~ 7 expected background events) 
with vanishing systematic uncertainty, which unfortunately is never the case in the environments at 
colliders. For the exclusion of a signal hypothesis, neglecting systematic uncertainties and performing a 
counting-only experiment with negligible background, a minimal signal yield of 3 expected events after 
application of the full selection chain is needed in order to exclude the signal+background hypothesis 
at the 95% confidence level, when using the Bayesian or "CLs" technique (see below). In this case 
the uncertainties on the signal rate and shape of the discriminating distributions will also influence the 
derived exclusion limits; for instance, a too optimistic assumption for the mass resolution for Higgs- 
boson candidates would yield a too stringent exclusion. 

Although the Higgs mechanism had been incorporated in the electroweak theory in the late 1960s, 
it lasted until 1989, when LEP started its operation, before significant parts of the mass range of the 
SM Higgs boson could be investigated^ The particle data group in 1988 [91] concluded on Higgs- 
boson searches before the start of LEP "In summary, the only cast-iron constraint on the Higgs mass 
is Mh > 14MeV. A combination of theoretical arguments and bounds from B, T, and K decays 
probably excludes the range below 4 GeV." Hence the mission of LEP was to probe the mass range of 
the SM Higgs boson starting from a zero mass to the highest possible mass values. For a discussion 
of Higgs-boson searches before the start of the LEP accelerator we, e.g., refer to Refs. [95T499] . In this 
article we review the most important results of the searches for the SM Higgs boson at the colliders 
LEP, Tevatron, and LHC. For details of the individual analyses, in particular identification criteria and 



6 A search carried out by the CDF experiment at Tevatron with up to 4.4 pb -1 at CM energy of 1.8 TeV excluded 
the mass ranges 2m p — 818 Me V and 846 MeV— 2mx at 90% CL [S3]. However, in this mass range the prediction of the 
branching ratios have large theoretical uncertainties, as discussed in Section [2~B1 
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systematic uncertainties, the original literature should be consulted. 
3.2 Statistical interpretation of search results 

The result of a search for the Higgs boson after applying the whole selection algorithms is, in the 
simplest case, given by the observed event yield, and the expected event yield for the background 
processes and the signal process, which depends on the hypothetical Higgs-boson mass. Such analyses 
are often called "counting experiments". Most of the searches performed during the LEP1 era were 
of this type. Very often the sensitivity of the search can be increased by using a final discriminating 
observable, which can be one-dimensional or, in rare cases, is chosen to be multi-dimensional. Here not 
only expected and observed event yields, but also the observed and expected shapes of the distributions 
in the final observable are used, when deriving the statistical interpretation of the results. This final 
discriminant is often given by the mass distribution of the Higgs-boson candidates, reconstructed from 
the decay products, or by the response of a multivariate event classifier (e.g. likelihood, artificial neural 
network, boosted decision tree). All searches during the LEP2 era, at Tevatron, and at the LHC 
use the distribution of a discriminating observable. Usually the expected event yields from signal 
and background processes as well as the shape of the distribution of the discriminating observable 
are affected by several systematic uncertainties, comprising, e.g., predictions of cross sections and 
acceptances due to kinematical requirements by theoretical calculations, knowledge of the detector 
performance in reconstruction and identification efficiencies, energy resolutions, and energy scales. Their 
influence on the signal and background expectations can be described by nuisance parameters 9, which 
can be common in signal and background. The expected distribution for a Higgs-boson signal with 
SM strength s(M H , 6) depends on the Higgs-boson mass M H and can either be given by a simple event 
count, a distribution of the final discriminant provided as a histogram, or an analytic function. In 
the following we will suppress the Mr dependence in the equations, but it should be kept in mind. 
The signal-strength parameter \i = cth/o"h sm allows to describe a signal strength different from the SM 
expectation, where /i = 1 holds. In the same way the background expectation is given by b{6) and the 
data by n. The statistical interpretation uses the frequentist concept of hypothesis testing for claiming 
evidence or discovery. When setting exclusion limits either the concept of hypothesis testing is used as 
well, or the closely related frequentist derivation of an upper limit on the signal yield or signal-strength 
parameter fi. For setting limits the Bayesian concept of providing a one-sided credibility interval or 
Bayesian upper limit is used as well. Below a partial overview of the main methods used at LEP, 
Tevatron, and the LHC is given. For more details we refer to the excellent review [100J and to the 
experiments' publications. A general introduction into the topics of frequentist hypothesis tests, and 
Bayesian and frequentist limit setting and their interpretation can be found, e.g., in Ref. [101J. 

In a hypothesis test the consistency of the observed data with a given hypothesis, which one tries to 
falsify (null hypothesis H ), is evaluated via the so-called p- value. The null hypothesis is the opposite 
of the statement one is aiming at, because hypotheses can only be rejected, but not be approved: 
When looking for evidence or discovery, it is given by the existence of only background processes 
called "background-only" hypothesis ("6" -only), whereas when excluding a signal, it is given by the 
"signal+background" hypothesis ("s + 6"). A test statistic t is constructed from the observed results 
(event yields and shapes of final discriminating observables) , which is used to evaluate the consistency 
with the null hypothesis. Given that the probability density function /(t(/z)|H ) for the test statistic 
t under the null hypothesis is known, and knowing that small (large) values of t indicate consistency 
(inconsistency) with the null hypothesis the p-value for an actual observed value t Q b s is defined by 



The p-value is the probability to observe a value of the test statistic t at least as large as observed in 
data assuming the null hypothesis is realized in nature. Equivalently it quantifies the probability to 
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wrongly reject the null hypothesis, if it is actually true. The p-value can be translated into a so-called 
Gaussian significance Z via Z = — p), where $ _1 is the inverse of the cumulative distribution for 

the standard Gaussian probability density function. For claiming discovery one requires the p- value to 
be less than or equal to 2.87 x 1CT 7 or equivalently the significance Z to be at least 5. For exclusion, 
when one wants to reject the signal+background hypothesis the p-value is usually required to be less 
than 0.05. The hypothesis is then called excluded with a confidence level CL = 1 — p. Alternatively one 
can determine the signal strength for which the p- value is exactly 0.05. This is equivalent to constructing 
a one-sided frequentist confidence interval for \i or an upper limit on /i called ^95 to a confidence level 
CL = 95%, which fulfills the condition 

0.05 = 1-CL=/ /(t(// 95 )|Ho)dt, i.e. 0.95 = CL = / /(t(/i 95 )|H ) dt. (29) 

All values of \i larger than /i 95 are then excluded, i.e. the SM Higgs-boson hypothesis is excluded if 
^95 < 1- The Bayesian technique to derive a one-sided credibility interval for a Bayesian upper limit 
on fj, at credibility level CL relies on the likelihood function £(data|/i), which gives the conditional 
probability to observe the data for a given value of /i. Using Bayes' theorem the posterior probability 
for the degree of belief in a particular value of \x in the light of the data is given by 

p(/i|data) = const, x £(data|//) 7r(/i), (30) 

where 7r(/x) is the a priori belief about the knowledge of /i and const, is an irrelevant constant to 
normalize p(/i|data). Very often and in all results discussed in this review, 7r(/z) is chosen to be uniform 
for non-negative and hence physically meaningful values of /i and zero otherwise. The Bayesian upper 
limit /igs to a credibility level CL = 95% is then obtained via 

Q g5 = a P fc|data)d/z = JT £(data|/x) d/i 
/-^ooP(/- t |data) d/i f£° £(data|/x) dfx 

The procedure described above is repeated for each hypothetical value of Mh in all approaches (fre- 
quentist hypothesis test, frequentist upper limit, and Bayesian upper limit). 

Let us first consider a counting experiment as was done during the LEP1 era and partially at 
the beginning of the LEP2 era. The discussion will be restricted to rejecting the signal+background 
hypothesis and hence to the derivation of exclusion limits, as during that era no significant excess was 
observed. The background b is assumed to be known without uncertainties and the expected signal 
yield is given by s. The upper limits on the signal yield is denoted by S95 = //95s. In the frequentist 
approach one has to specify the test statistic, which is chosen to be the number of observed events n, 
which in this case does not depend on the signal hypothesis. The probability density function for n 
under the signal+background hypothesis is simply the Poisson probability f(n\pLS + b) = e~^ J,s+h ' (p,s + 
b) n /n\. The likelihood function needed in the Bayesian approach is also given by the Poisson probability 
£(data|yus + b) = e~^ s+b \fis + b) n /n\, where n plays the role of "data". A low number of observed 
events indicates an inconsistency with the signal+background hypothesis. Hence the frequentist upper 
S95 limit is derived from 

^■obs 

0.05 = J2 e- {s95+b) (s 95 + b) n /n\ . (32) 

n=0 

The analytic solution of this equation for s g5 gives 

s 95 = ^(0-95; 2(n obs + 1)) - 6, (33) 

where F~2 denotes the inverse of the cumulative x 2 probability density distribution evaluated for 
2(n Q bs + 1) degrees of freedom. The Bayesian upper limit has to be derived from [102] 

n n , _ Sn= b o e-^ +b \s 9 , + b) n /n\ 



23 



For zero background {b = 0) the two approaches give identical results. For n bs = (1) observed events 
the S95 is then given by 3.0 (4.7). These are the numbers that are used in the derivation of the final 
LEP1 exclusion limits, where all observed events were considered as candidate events or equivalently 
using b = in the above equations for S95. A Higgs-boson mass hypothesis is excluded if the expected 
signal event yield for the value of Mh tested is larger than 3.0 (4.7) given (1) events are observed 
in data. The expected signal event yields are decreased by their systematic uncertainty before the 
determination of S95 and its comparison with [is. 

For non-vanishing background the frequentist approach due to the —b term can give arbitrarily 
small values S95, e.g. S95 = for 6 = 3 and n = 0, if the observed event yield is significantly lower 
than the one expected from background processes. In the Bayesian approach S95 can never get smaller 
than 3. Even if completely correct from a purely frequentist point of view, people felt uneasy with 
quoting a very small S95 in a situation where the average experiment has no sensitivity to exclude 
such a small signal yield, but excludes it due to a downward fluctuation of the observed event yield 
with respect to the background-only hypothesis. Such situations are avoided in the Bayesian approach. 
An alternative derivation of the Bayesian formula in a pseudo-frequentist ansatz and its interpretation 
are given in Ref. |103j . Formally the right-hand side of the Bayesian formula can be rewritten in the 
language of frequentist hypothesis tests as p s+ b/ (1 — Pb), where p s+ b and p^ denote the p- values under the 
signal+background and background-only hypotheses, respectively, keeping in mind that small (large) 
observed event yields mean inconsistency with the signal+background (background-only) hypothesis. 

At LEP2 and later, multiple channels with different signal-to-background ratios were combined, 
and the distributions of the final discriminants were also used. At the beginning of the LEP2 era 
different techniques |104j were used for the construction of a test statistic that takes into account the 
different signal-to-background ratios in the various search channels and bins of the histograms of the 
final discriminants. The likelihood functions £(data|/xs + b) and £(data|6) are given by the product 
of Poisson probabilities running over all channels and bins. Neglecting systematic uncertainties, the 
optimal test statistic is given, according to the Neyman-Pearson lemma |105j . by the ratio of the 
two likelihoods to observe the data under the signal+background hypothesis and background-only 
hypothesis, its reciprocal or any monotonic function of it. This is the statistic 



or its transform —2 InQ used at LEP. For a single-channel/single-bin counting experiment the Neyman- 
Pearson test statistic derived from the ratio of two Poisson probabilities is (setting jj, — 1) 



Hence as Q is monotonic in n, n is the optimal test statistic for a simple counting experiment in one 
channel. Each observed event enters the calculation of —2 InQ with a weight of ln(l + s/b) which takes 
into account the local s/b ratio. Large values of —2 InQ are obtained for "b-only"-like data and small 
values for "s+b"-like data. The median of the —2 In Q distribution is positive under the background-only 
hypothesis, whereas the median is negative under the signal+background hypothesis. The probability 
density functions for —2 In Q under the signal+background hypothesis and background-only hypothesis 
are obtained from simulation of toy MC experiments. 

Translating the Bayesian formula in the interpretation of Ref. |103] to the test statistic —2 In Q 
yields the so-called "CLs" technique |106l H07j. which avoids exclusion without sensitivity. The p- values 
for the LEP test statistic —2 In Q are calculated according to: 



Q = £(data|/is + 6)/£(data|6), 



(35) 



Q = e~ s (l + s/b) n , 



i.e. 



21nQ = 2s-2nln(l + s/6). 



(36) 




(37) 



and 




(38) 
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The "CLs" value is then calculated as 



CLs=p s+ b/(l-Pb). (39) 

A Higgs-boson mass hypothesis is called excluded at least to 95% CL if CLs < 0.05. The probability 
to reject the signal+background hypothesis, if it is true, is less than 5%, as the "ad hoc" correction in 
the denominator destroys the strict frequentist interpretation. Even if not proven in a mathematically 
rigorous way, good agreement has been found between the Bayesian approach and "CLg" technique in 
practice at LEP, Tevatron, and the LHC. Systematic uncertainties (denoted by 0) were incorporated 
at LEP in the hybrid frequentist-Bayesian approach described in Ref. [108J, which results in modified 
probability density functions for the unchanged test statistic —2 In Q, 

f(-2lnQ\s + b) = J f{-2\nQ\s + b,6)7i{6)de and f{-2 In Q\b) = J f{-2 \nQ\b, 0)tt(0) d0, (40) 

which is a so-called marginalization with respect to the parameters with a priori probabilities 7r(0). 
The explicit choice of 7r(0) is specific to the sources of the systematic uncertainties and thus depends 
on the analysis. 

At Tevatron the Neyman-Pearson test statistic was modified to include the parameters 6, describ- 
ing the systematic uncertainties, in the likelihood, which yields a better discrimination between sig- 
nal+background and background-only hypotheses in the presence of systematic uncertainties compared 
to LEP. The likelihood function now reads £(data|/is(0) + b(9))n(9), with a priori probabilities 7r(0) for 
the nuisance parameters 0, and the test statistics LLR, in the Tevatron convention, is given by ratio of 
two profiled likelihoods 

LLR = -2 1 nf £ ' data '^' J > + ^^). (41) 
\ £(data|6(0))7r(0) J 

The and are, respectively, the conditional maximum-likelihood estimators of two fits performed (i) 
under the signal+background hypothesis in the numerator and (ii) under the background-only hypoth- 
esis in the denominator. The best fit values and for 8 are usually different under the two hypotheses. 
The determination of the probability density functions for LLR and of the "CLs" value proceeds as 
was done at LEP. The primary choice for deriving upper limits on /i, which is called R in the Tevatron 
analyses, in the latest Tevatron combinations is the use of the Bayesian method with the posterior 
probabilities for /i given by 

p(/x|data) = const, x J £(data|^s(0) + b(9))ir(9)ir(Li) d0. (42) 

Again a uniform prior for positive p is used and a marginalization with respect to the is performed. 
For deriving p-values under the background-only hypothesis and calculating significances, a fit of 
£(data|/is(0) + b(9))ir(6) to the data is performed, and the best fit value p for the signal-strength 
parameter p is used as the test statistic. The probability density functions for p are derived again from 
toy MC experiments. 

At the LHC two different test statistics are used for the two cases, setting limits and claiming evidence 
for discovery |109j . The likelihood function is the same as at Tevatron, £(data|/is(0) + b{9)) p{9\9) , but 
the a priori probabilities for the nuisance parameters 7r(0) are reinterpreted or replaced by conditional 
probabilities p{9\9) to observe a value in an auxilary measurement given the value 0. For setting 
limits, i.e. testing a signal strength p, the test statistic is given by 

2 ^( c(d a t a \ps(§) + b(§))p(9\§) \ 
\£(data|/xs(0) + 6(0))p(0|0) ) ' 
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In the numerator the likelihood under the signal+background hypothesis with fixed signal strength /a, 
to be tested, is evaluated, while in the denominator (in contrast to Tevatron) the likelihood is evaluated 
after fitting also for the best signal strength consistent with the data ft. The range for fx is restricted 
to < p. < fi, as fx < is unphysical and as a one-sided hypothesis test is performed, where fi > ll 
is considered as being consistent with the signal+background hypothesis (see Ref. |109] for details). 
The value of the test statistic g M is restricted to the range between zero and infinity. In order to be 
able to define a one-sided hypothesis test, the value of g M is set to zero for fi values larger than ll. 
Background-like data yield large values of g M , whereas "s + 6" -like data accumulate at low values of 
g M . The advantage of not fixing ll to in the denominator, as done at LEP and Tevatron, is that 
the probability density function for in the limit of large event yields are known analytically due to 
Wilks' [110] and Wald's [1 llj theorems (see Ref. |109j for the approximative analytic expressions). The 
determination of frequentist upper limits proceeds as at LEP and Tevatron via the "CLs" technique. 

For testing the background-only hypothesis, when looking for discovery or evidence at the LHC, the 
test statistic is calculated as a different ratio of profiled likelihoods: 

* = -*J Whm\. (44, 
\£(d&t&\fis{9) + b(9))p(9\9)J 

Now ll is fixed in the numerator under the null hypothesis (background-only) and is determined via 
a fit to the data fi in the denominator. The range for fi is restricted to > 0, as ll < is unphysical 
and as a one-sided hypothesis test is performed, where fi < is considered as being consistent with 
the background-only hypothesis (see Ref. [109J for details). The values of go are restricted to the range 
between zero and infinity. In order to be able to define a one-sided hypothesis test, the value of q is set 
to zero for fi < 0. Now "6-only"-like data accumulate at low values of g , whereas "s + 6"-like data yield 
large values of go- Again the test statistic go has the advantage that its probability density function is 
know analytically due to Wilks' and Wald's theorems, i.e. 

/(goi&) = ^%o) + ^(9o;l), (45) 

by a combination of ^-function and |x 2 -P roD ability density function with 1 degree of freedom (see 
Ref. |1U9] for details). The |<5-function is due to the fact that go is set to for fi < 0. Hence any 
observed go ^ yields a maximal p- value of 0.5. 

So far in the likelihood functions that enter the construction of test statistics used for testing the 
background-only hypothesis (— 21nQ at LEP, fi at Tevatron, g at the LHC) the Higgs-boson mass 
hypothesis Mh is fixed under the alternative signal+background hypothesis. The corresponding p- 
values are hence called local p- values. It gives the probability to observe such a fluctuation for a fixed 
Mh value. The so-called global p-value denotes the probability to find such an excess anywhere in 
the mass spectrum or more precisely for arbitrary Higgs-boson mass hypothesis under the alternative 
hypothesis within a certain Mh range. The relation between the two p-values is given by p g iobai = /trial 
Piocai- The trial factor / tria i can be interpreted as the effective number of independent searches performed. 
This enhancement in the p- value is commonly called the "look-elsewhere effect" . A technique to derive 
approximate global p- values and trial factors used at the LHC can be found in Ref. |112j . and in 
Ref. [113j for Tevatron. 

Often expected limits on the signal strength, expected p-values, and significances are quoted. These 
are the median values for those quantities derived under the alternative hypothesis: Expected limit is 
the median limit under the background-only hypothesis, expected significance is the median significance 
under the signal+background hypothesis. 

A subtle difference should be noted when generating toy MC experiments at Tevatron and the LHC. 
The central values of the nuisance parameters, which are used as input to the MC experiments, are 
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fixed to the a priori best values at Tevatron and are extracted from a fit to the data at the LHC (see 
Ref. [333 for details). 

In the following, all limits are given to 95% CL, if not stated otherwise. 

4 Higgs production at e + e _ colliders 

At e + e~ colliders the importance of the various Higgs production mechanisms strongly depends on the 
centre-of-mass (CM) energy of the collider and to some extent also on the Higgs-boson mass. Figure [7] 
illustrates the production channels by some representative lowest-order diagrams. Up to LEP energies, 
y/s < 209 GeV, ZH production via "Higgs-strahlung" with a variety of Higgs and Z-boson decays was 
the main search channel. Some events from H7 production were expected in spite of a very small 
cross section; owing to the extremely suppressed electron Yukawa coupling the leading order of the 
channel is loop-induced by massive-fermion and W-boson loops. At future e + e~ linear colliders, such 
as the ILC or CLIC, two additional types of production channels become relevant: the fusion of W 
or Z bosons and the associated production with tt pairs. For instance, the cross section of W-boson 
fusion at an ILC with a CM energy yfs = 500 GeV exceeds the one of Higgs-strahlung off Z bosons 
for Higgs-boson masses larger than ~ 160 GeV. While the HZ production cross section falls off oc 1/s, 
like any typical s-channel process, the t-channel-like W/Z-fusion processes receive contributions that 
grow like ln(s / M^) / for large energies. Above the ttH production threshold, y/s > 2m t + Mh, 
Higgs production in association with tt pairs offers the possibility to directly measure the top-quark 
Yukawa coupling with good precision. Figure M (r.h.s.) presents an overview over the various cross 
sections as function of the collider energy y/s, covering the range from LEP energies to the TeV region 
of potential future colliders. The theoretical preparation of the LEP analysis is described in detail in 
Refs. |115[I118] , and a more detailed survey of both the physics potential of linear colliders and features 
of the predictions can be found in Refs. [0*9 11119] . 

Figure M is based on lowest-order predictions, but state-of-the-art calculations include higher-order 
improvements as well. Collinear photonic initial-state radiation (ISR) generically delivers corrections at 
the level of 10% or more, which are particularly pronounced in region where cross sections or distribu- 
tions show strong variations, as it is, e.g., the case near kinematical thresholds. The by far dominating 
part of ISR is of universal origin and consists of logarithmically enhanced terms a n ha m (m e /Q) with 
(m < n), where Q is a typical scale of the process. This important effect can be taken into account by 
structure functions, which are known to quite high order (see, e.g., Ref. |120] and references therein). 
The full inclusion of NLO corrections in predictions, however, require dedicated calculations: 

• Higgs-boson-photon production e + e~ — > Hy [121U122j 

The production of a Higgs boson in association with a photon mainly proceeds via top-quark 
and W-boson loops, where a significant destructive interference occurs. At LEP energies, H7 



Figure 7: Representative lowest-order diagrams for the main SM Higgs-boson production channels at 
e + e~ colliders: (a) Higgs-strahlung, (b) vector-boson fusion, (c) top-quark associated production, (d) 
Hy production. 
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Figure 8: Left: cross section for e + e~ — > Hvv for a typical LEP2 energy and its breakup into Higgs- 
strahlung, W-fusion, and interference contribution (taken from Refs. [115(1116] ); right: cross sections for 
the various SM single-Higgs-boson and Higgs-boson pair production channels for a future e + e~ linear 
collider (taken from Ref. |117j ). 

production only contributed in the search for a light Higgs boson with a mass well below the Z- 
boson mass, i.e. H7 production was mainly treated as potential decay mode of the Z boson [121] . 
At a high- luminosity e + e~ linear collider H7 should, however, also be observable at energies far 
above the Z resonance |122] . where e + e~ — > H7 is a true 2 — > 2 scattering process. Since this 
rare process is already loop induced, full NLO predictions would require a two-loop calculation 
with multiple scales that is technically still out of reach with present calculational techniques. On 
the other hand, taking into account the universal ISR effects most likely will also be sufficient to 
match the expected experimental precision at a future collider. 

• Higgs- strahlung e+e" HZ [5T1IT21M21)] 

The actual Higgs-strahlung signal that was searched for at LEP was mainly based on the H — > 
bb decay and various Z-boson decays Z — > ff. Higher-order effects from ISR lead to relevant 
corrections of 10—20% at LEP energies and have been taken into account in the LEP analyses. 
The remaining genuine weak corrections, which comprise weak-boson exchange etc., have been 
first calculated in the approximation of on-shell Z and Higgs bosons |124] and amount to 1—2% 
only at LEP energies. Later the EW corrections have also been evaluated including the leading 
off-shell effects of Z — > fi + fi~ and H — > bb |125j . which will be relevant at possible precision studies 
at a future linear collider. The EW corrections become more and more relevant with increasing 
collider energy. 

• W fusion e + e"(WW) -> Hw [TT6lfT2^[TMfl3T] 

Actually the physically observable process involves a sum over all three neutrino species in the 
final state. While W fusion is only possible for v e v e pairs in the final state, all three vv channels 
receive contributions from s-channel HZ production with Z — > vv decays; for B.v e v e t- and s- 
channel interfere. The l.h.s. of Figure [8] illustrates the interplay of HZ(— > vv) production and W 
fusion for a typical LEP2 energy near the HZ threshold region where W fusion becomes significant. 
For -y/s <• Mh + Mz the HZ production process can only proceed via off-shell Z bosons (the Higgs 
boson is too narrow to get off shell for the Mh values relevant at LEP). Sufficiently above the HZ 
threshold (y/s > M H + M z + 100 GeV for M H ~ 100-200 GeV), which is an interesting region 
for future linear colliders, W fusion dominates (see r.h.s. of Figure [8]). The NLO corrections are 
of pure EW origin and known [1301 1131] . Owing to the dominating t-channel kinematics, which 
prefers forward/backward-produced (anti)neutrinos and low-virtuality intermediate W bosons, 
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Figure 9: Lowest-order diagrams for Higgs-boson pair production at e + e colliders: (a) Higgs-strahlung, 
(b) W/Z-boson fusion. 

the corrections to W fusion can be described within few percent by an approximation based on 
universal ISR (with appropriate scale Q) and a simple constant correction factor from heavy top 
loops [131] . 

• Z fusion e+e"(ZZ) He+e" [T27HT29l[T32] 

The Z-fusion cross section is suppressed with respect to the one of W fusion by roughly a factor 
of 10 as a mere consequence of the different size of the charged- and neutral-current couplings of 
the electron. However, Z fusion bears the advantage that it also would be observable at the e~e~ 
variant of a linear collider. Similar to HVz/ production, He + e~ production receives contributions 
from s-channel HZ production diagrams and t-channel Z fusion graphs, where the logarithmic 
rise of the latter dominates at high energies. The EW corrections [132] again receive only small 
contributions < 5% beyond ISR. 

• tt associated production e + e~ — > ttH [133H140] 

A precise measurement of the top-quark Yukawa coupling represents the main motivation to an- 
alyze ttH production in the clean environment of an e + e~ collider [134] . Assuming, for instance, 
M H = 120 GeV the Yukawa coupling could be measured with an accuracy of ~ 5% at a linear col- 
lider operating at the energy of 800 GeV after collecting the integrated luminosity of 1 ab _1 [135J; 
in combination with other measurements even a higher accuracy is possible [136] . Owing to the tt 
pair in the final state, the ttH cross section receives both QCD [137] and EW corrections [138] at 
the NLO level, which are completely known. Sufficiently above the tt threshold the generic size 
of QCD, photonic, and remaining electroweak corrections is 5—10% on the cross section, but the 
effect on kinematical distributions can be much larger. In particular, QCD correction in the tt 
threshold region deserves particular attention, since it is dominated by the Coulomb singularity. 
A careful treatment of these long-distance effects within non-relativistic QCD [139] reveals cor- 
rections of ~ 70% relative to the NLO cross section for slowly moving tt pairs. For a Higgs boson 
with M H = 120 GeV at a 500 GeV collider this pushes the NLO K factor of 1.7 further up to 2.4. 

On the other hand, a realistic, precise description of the process has to include the top-quark 
and Higgs-boson decays, leading to a 10-particle final state already at LO, assuming the Higgs- 
boson decay into bb pairs. Corresponding multi-particle Monte Carlo simulations are described 

mRef. ma. 

A future e + e~ linear collider could also allow for a determination of the triple Higgs-boson self- 
coupling via Higgs-boson pair production [117jll4lT - |143] if the collider energy is high enough. As 
illustrated in Figure [91 the promising production channels are Higgs-strahlung and W fusion, where 
diagrams with triple Higgs couplings compete with other pair production diagrams. Figure [8] includes 
the corresponding Higgs pair production cross sections for a linear collider up to energies in the few-TeV 
range. Similar to single-Higgs-boson production, the Higgs-strahlung mechanism delivers the largest 
production cross section for Higgs-boson pairs for energies directly above the HHZ threshold, but this 
cross section decreases oc 1/s. W fusion again shows a logarithmic rise with increasing energy and takes 
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over the leading role for energies y/s > 1 TeV for Mh = 120 GeV. It should be noted, however, that for 
a low Higgs-boson mass the cross sections are < 0.2 ftT 1 for y/s < 1 TeV and get even more suppressed 
for larger Mh. Thus, a determination of the triple-Higgs coupling to ~ 20% for Mh = 120— 180 GeV 
requires a high-luminosity linear collider (L ~ lab -1 ) operating at TeV energies [142J. The EW 
corrections both to HHZ production [143 J and to W fusion to Higgs-boson pairs |117j are known. The 
results resemble the qualitative features of their single-Higgs counterparts, but the size of the effects get 
somewhat larger. In addition to the obligatory large ISR effects, the genuine weak corrections are at 
the generic level of 5—10%. For s-channel HHZ production the weak corrections show the tendency to 
grow negative with increasing energy, while the weak corrections to t-channel W fusion hardly depend 
on the energy sufficiently above the HH threshold. 

A direct experimental investigation of the quartic Higgs self-coupling via triple Higgs production 
processes seems out of reach, since the cross sections are too much suppressed and are not very sensitive 
to the H 4 coupling. 

5 Searches at LEP 
5.1 The LEP1 era 

From autumn 1989 to summer 1995 the LEP accelerator was operated at CM energies in the vicinity 
of the mass of the Z boson, and in total 15.5 million hadronic Z-boson decays Z — > qq and 1.7 million 
leptonic decays Z — > £ + £~ were collected by the four experiments ALEPH, DELPHI, L3, and OPAL. 
At the end of LEP1 data taking a peak luminosity of 2 x 10 31 cm _2 s _1 was reached, corresponding to 
approximately 1000 Z bosons being recorded per hour by each experiment. Even for very low Higgs- 
boson masses very efficient triggering of the final states of interest was possible and did not impose a 
severe problem for the analysis. 

As detailed in Section HJ the dominant production mode for the Higgs boson is the rare decay 
Z — > Z*H. The total expected event rate for the decay Z — > Z*H — > /j/jH for one million of hadronic 
Z decays is approximately 75 (12,0.4) for a Higgs-boson mass of 10 (30,60) GeV (see Figure [TUl left). 
For a Higgs-boson mass of roughly 60 GeV the event rates in Z — > Z*H — > /j/jH and Z — > 7H are 
of comparable size (see Figure [TDj left). However, considering also the decay of the Z* into electrons 
and neutrinos increases the rate in accessible final states arising from Z — > Z*H by a factor of eight. 
In addition to the smaller rate for Z — > 7H the background from ffj was also more severe for this 
final state. Hence no dedicated search for the decay Z — > 7H was performed, as the contribution to 
the sensitivity is negligible given the integrated luminosity collected at LEP1 of up to approximately 
160 pb _1 . At LEP1 the Higgs-boson mass range from to approximately 65 GeV could be probedQ The 
a priori signal-to-background ratio corresponds to 10~ 2 for low masses and 10~ 6 at the edge of the mass 
sensitivity (see also Figure [6]). The kinematics in Higgs-boson production and decay is characterized by 
the average values of the Higgs-boson momentum and production angle, which are shown in Figure [10] 
(middle and right) as function of the Higgs-boson mass. 

Three Higgs-boson mass ranges can be roughly distinguished: (a) < Mh < 2771^, (b) 2m M < Mh < 
3 GeV, and (c) Mh > 3 GeV. These mass ranges differ in the dominant Higgs-boson branching ratios 
and the robustness of their theoretical predictions, the kinematics and topology, the composition of the 
dominant backgrounds, and the search strategies applied. 

• For 2m M < Mh < 3 GeV a perturbative description of the Higgs-boson decay, the prediction of the 
branching ratios, and hence the estimation of the detection efficiencies have large uncertainties. 

7 In parallel to the start of LEP a search for light Higgs bosons, performed in an electron beam-dump experiment, was 
published, which unambiguously excludes the mass range 1.2— 52MeV [144] . 
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Figure 10: Relative decay rates for Z — >■ Z*H — > /i/nH and Z — > 7H and number of expected events per 
10 6 hadronic Z decays [145J (left), average momentum of the Higgs boson (middle), and average opening 
angle of the Z* and H-boson decay products [98J (right) as a function of the hypothetical Higgs-boson 
mass. 



The Higgs boson decays dominantly into several hadrons (e.g. a pair of pions, kaons, etc.) and for 
2m^ < Mr < 2777,7,- into a pair of muons. Different approaches were used by the four experiments 
to cover this mass range. ALEPH |146j and DELPHI [147] extended their standard searches 
in the lepton Z* — > £ + £~ and missing energy Z* — > vv final states (see below) to the mass 
range below 3 GeV, allowing only two tracks being assigned to the Higgs-boson candidate. They 
relied on some modelling of Higgs-boson decays in this non-perturbative regime. In addition, 
for mass values below 2 GeV the decay Z* — > qq was considered by asking for three- or four- 
jet events and assigning the jet with the lowest track multiplicity to the Higgs-boson candidate, 
which is only allowed to have two tracks. L3 [148] only utilized the decay Z* — > £ + £~ with 
two energetic acoplanar leptons and considered Higgs-boson decays to a least one track, covering 
the decay into muons and hadrons. OPAL [149J performed a decay-mode-independent search 
by two complementary analyses: (a) acoplanar leptons with vetoing on H — > "electromagnetic" 
e + e _ /77/7r 7r° decays, which is also sensitive to a long-lived Higgs boson decaying outside the 
detector, and (b) a dedicated search H — > "electromagnetic" decays accompanied by missing 
energy from Z* — > vv. The OPAL approach does not rely on the prediction of branching ratios 
and details of modelling the hadronic Higgs-boson decay in the non-perturbative regime, and 
hence is almost model independent. 



• For Mh < 2777^ the dominant decay is to a pair of electrons, and for Mh < 2777 e only the loop- 
induced decay to a pair of photons is possible. In this mass range the Higgs boson gets long-lived 
(average decay length of 1 (100) meter(s) for Mh = 100 (10) Me V), and a significant fraction 
of Higgs bosons decay outside the detector, leaving an invisible Higgs-boson decay as signature. 
ALEPH [M[T50] and DELPHI [151] extended the above searches in the Z* ->• qq and Z* vv 
final states and performed a dedicated analysis, which looked for a V° signature (i.e. a displaced 
vertex with two tracks) accompanied by any decay mode of the Z* boson. The di-track invariant 
mass of the V° candidate should not be consistent with that of a long-lived hadron (K®, A, etc.), 
and the tracks should be not identified as hadrons. L3 [152] used the acoplanar lepton search from 
the 2777^ < Mh < 3 GeV range allowing also for H — > e + e~ decays. OPAL [153J looked for Z* — > vv 
and H — > e + e~ /jj decays with at least one energy deposit in the electromagnetic calorimeter, 
to which tracks can be assigned, and not much additional detector activity. All experiments 
complemented this analysis by a search for invisible Higgs-boson decays in the acoplanar lepton 
topology with no other significant activity in the detector. A combination of the "invisible Higgs- 
boson" search and the other ones discussed above exclude all mass values down to 0. 
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• For Mh > 3 GeV the perturbative description of the Higgs-boson decay is more reliable (apart 
from the mass values corresponding to bb and cc bound states). For low Higgs-boson masses 
(Mh < 10— 15 GeV) the Higgs-boson decay products build a mono-jet topology due to the large 
Lorentz boost. For larger Higgs-boson masses the topology tends to two separated jets, which 
have an opening angle less than 180 degree in space ( "acollinearity" ) and also in the projection 
into the plane perpendicular to the beam axes ( "acoplanarity" ) . The final state with largest 
sensitivity is the decay Z — > Z*H — > vvH, where for Mh larger than 15 GeV the decay into a pair 
of b-quarks is dominant. A production rate that is smaller by a factor three is provided by the 
decays Z — > Z*H — > yU + yU _ H/e + e _ H. This ratio of three is roughly retained at the end of the event 
selections. Together the two final states cover approximately 25% of the Higgs-boson production 
rate. Final states involving tau leptons, Z*H — > r + r~qq, Z*H — > qqr + r~ , Z*H — > £ + £~t + t~, 
and Z*H — > vvt + t~, have either a small production rate or a quite low selection efficiency The 
expected signal event yields after the full selection are approximately one order of magnitude 
smaller than that expected in the Z*H — > vvW. search. Hence final states with tau leptons only 
contribute marginally to the combined sensitivity. All collaborations investigated final states with 
tau leptons [150| I154H156] at intermediate stages of the Higgs-boson search at LEP1, but not 
included them in the analysis providing the final mass limit. Fully hadronic final states from 
Z — > Z*H are not considered due to the overwhelming background. 

The Z Z*H — > vvH topology is characterized by missing energy, missing (transverse) momen- 
tum, and the detector activity from the Higgs-boson decay. With increasing Higgs-boson mass 
the particle multiplicity from the Higgs-boson decay is increased. Due to the smaller Lorentz 
boost the topology changes from a mono-jet to an acollinear and acoplanar di-jet final state. 
At the same time the missing energy and missing mass decrease, as the Z* boson has to be 
more off shell. Almost no irreducible backgrounds exist, and the background contributions are 
mostly of instrumental nature due to mismeasurements, wrong identification of final-state parti- 
cles, and incomplete detector coverage. The dominant background for large Higgs-boson masses 
is due to Z — > qq decays, where missing energy arises from neutrinos produced in heavy-flavour 
decays and mismeasurements, and from four-fermion production via Z* — > Z*7* — > vvqq and 
e + e~ — > e + e~7*7* — > e + e~qq. The four-fermion processes are characterized by a low mass of 
the hadronic system and hence mostly relevant for small Higgs-boson masses, and events of the 
e + e~qq channel typically show small (missing) transverse momentum. The mass resolution for 
the hadronic system is of order of 10%. DELPHI also applied b-tagging for the final analysis 
exploiting the large branching ratio H — > bb. 

The Z — > Z*H — > £ + £~H topology is characterized by two isolated leptons, whose acollinearity and 
acoplanarity increase with increasing Higgs-boson mass, and a recoiling hadronic system. The 
mass of the Higgs-boson candidate can be reconstructed from the four-momenta of the leptons 
using energy-momentum conservation ("recoil mass technique") with a precision of 1—2%. The 
dominant background from four-fermion production has either a small mass of the leptonic or 
of the hadronic system. The Z — > qq background with leptonic heavy-flavour decays possesses a 
large mass of the hadronic system, and the leptons do not originate from the primary collision 
vertex. For the final analysis performed at LEP1 most experiments exploited b-tagging to further 
discriminate the signal especially from four-fermion background processes. 

The event selections were mostly designed in such a way that the expected background was at 
the level of one event or less at the edge of the Higgs-boson mass sensitivity for a given data 
set. With increased integrated luminosities this required more sophisticated and refined analyses 
and tighter identification and isolation requirements to exploit the above-mentioned differences 
between signal and background topologies and their dependence on the Higgs-boson mass. Hence 
the selection efficiency for signal events at the sensitivity edge dropped from roughly 60% or higher 
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Table 2: Excluded Higgs-boson mass ranges from LEP1 searches: experiment, excluded mass range, 
number (in thousands) of recorded hadronic Z decays, integrated luminosity, years of data taking and 
references. The entry "-" indicates that the information is not available. 



Experiment 


excluded -Mh range 


10 3 Z -> qq 


int. lumi. [ pb ] 


data from year 


Ref. 


ALEPH 


032-15 GeV 


12 


0.5 


89 


[146J 




11-24 GeV 


25 


1.2 


89 


[161] 




0-57 MeV 


23 


1.2 


89 


[162] 




< 41.6 GeV 


100 




89-90 


[150] 




< 48 GeV 


185 




89-90 


[163] 




< 58.4 GeV 


1233 




89-92 


[164] 




< 63.9 GeV 


4500 




89-95 


11571 


DELPHI 


2m„-14GeV 

P 


11 


0.5 


89 


|147| 




0-2m„ 

P 


119 




90 


[151] 




< 38 GeV 


119 




90 


[154] 




< 55.7 GeV 


1000 


34.6 


91-92 


[158J 


L3 


2-32 GeV 


50 




90 


[148] 




0-2 GeV 


70 




90 


[152] 




< 41.8 GeV 


111 


5.3 


90 


[165] 




< 52 GeV 


408 


17.5 


90-91 


[155] 




< 57.7 GeV 


1062 


39.0 


90-92 


|166j 




< 60.2 GeV 


3000 


114 


91-94 




[1591 




OPAL 


3—19 3 GeV 




0.8 


89 


[TB7] 




0-2m u 


25 


1.2 


89 


[153] 




3-25.3 GeV 


25 


1.2 


89 


[168] 




3-44 GeV 


170 


8 


89-90 


[156] 




0-11.3 GeV 




6.8 


90 


[149] 




< 56.9 GeV 


1900 


78 


90-93 


[169] 




< 59.6 GeV 


5000 


160 


89-95 


[160] 



in the first searches to a level of 30% in the final LEP1 searches. Only a few events survived all 
selection requirements in agreement with expectations from background processes. 

All observed events are considered as Higgs-boson candidates and the expected signal yields are 
reduced by their systematic uncertainty. In summary, the searches performed by the four experiments 
allowed to exclude all Higgs-boson mass hypotheses below 63.9 GeV (ALEPH |157j ). 55.7 GeV (DELPHI 
[158]), 60.2 GeV (L3 [159]), and 59.6 GeV (OPAL [IB0] ). as shown in Figure [LT] and in Tabled The 
relative contribution to the sensitivity of the missing energy and leptonic final states can be deduced 
from Figure [TTJ In ALEPH and L3 the mass resolution of candidates close to the edge of the excluded 
mass range is taken into account following the technique described in Ref. |170j . whereas OPAL did 
not apply such a technique. No official combination of the findings of the four LEP experiments was 
performed. Owing to the exponential decrease of the Higgs-strahlung cross section with increasing 
Higgs-boson mass, the excluded mass range only extended slowly with collected integrated luminosity 
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Figure 11: Excluded and expected event yields for Higgs-boson production at LEP1 by ALEPH (top- 
left) [157], DELPHI (top-right) [[58], L3 (bottom-left) [159], and OPAL (bottom-right) [T60] . 



(see Table E]). Hence an extension of data taking at LEP1 would not have significantly increased the 
mass reach of the searches, which for the first time excluded a significant fraction of the allowed Higgs- 
boson mass range in an unambiguous way. 



5.2 The LEP2 era 

Starting in summer 1995 at 130/136 GeV, the CM energy of the LEP collider was continuously increased 
over the next years up to 209 GeV in 2000. The data taken at 130/136 GeV were not used for SM Higgs- 
boson searches. A summary of the data sets recorded at different CM energies used for SM Higgs-boson 
searches during the LEP2 programme is given in Table |3j Approximately 0.5 (2.5) fb _1 at CM energies 
in excess of 206 (189) GeV were collected in total by the four experiments ALEPH, DELPHI, L3, and 
OPAL. At LEP2 CM energies the dominant production process is again Higgs-strahlung. But in contrast 
to LEP1 the Z boson in the final state is dominantly produced on its mass shell. The difference in the 
M H dependence of e+e~ Z ->• Z*H at LEP1 and e+e" -> Z* ZH at LEP2 is shown in Figured! 
(left). Additional contributions from the vector-boson-fusion processes e + e~ — > e + e~Z*Z* — > e + e~H and 
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Table 3: CM energies and integrated luminosities collected in each experiment at LEP2. Numbers in 
brackets give the luminosity-weighted CM energy of the corresponding year. The luminosity ranges are 
caused by the different data-taking efficiencies in the four experiments. 



year 


1996 


1997 


1998 


1999 


2000 


CM energy [GeV] 
int. luminosity [pb _1 ] 


161/170/172 
10-11/1/9-10 


183 
54-57 


189 
158-176 


192-202(198) 
228-237 


202-209(206) 
217-224 




Figure 12: Comparison of Higgs-strahlung production cross sections at LEP1 and at the beginning of 
LEP2 (left from Ref. [171] ). Number of expected Higgs-boson events per experiment produced in 2000 
data taking from Higgs-strahlung (full line) and from vector-boson fusion and interference (dashed line) 
(right from Ref. [172] ). 

e + e~ — > uuW*W* — > vvH and their interference with the respective final states from Higgs-strahlung 
are at the level of maximally 10—20% (see Figured]) at the sensitivity edge, which is approximately 
given by — M% + few GeV. The contribution to the total event yield (per experiment) for the 
production of a hypothetical Higgs boson in 2000 data is shown in Figure [12] (right). No dedicated 
searches for production in vector-boson fusion were performed. Signal cross sections were calculated 
and signal event were generated with HZHA [173J. 

In contrast to LEP1, now the hadronic decays of the Z boson can be utilized in the search, as the 
a priori signal-to-background ratio is larger by two to three orders of magnitude compared to LEP1 
(see l.h.s. of Figure [TO]) at the sensitivity edge. Furthermore the mass constraint of two jets to stem 
from the decay Z — > qq helps to discriminate the signal process from background processes. This new 
channel delivers the largest contribution to the sensitvity at LEP2. At LEP1 the background was 
mostly reducible, at least in principle, and caused by mismeasurements and hence difficult to model 
in all details. After crossing the thresholds for WW and ZZ production at the LEP2 CM energies of 
161 GeV and 183 GeV, respectively, the dominant backgrounds are due to di-boson production and to 
a large extent irreducible. These could reliably be predicted via simulation. 

Searches were performed in the "four-jet channel" ZH — > qqbb (ZH — > qqqq), "missing-energy chan- 
nel" ZH — > vvbb (ZH — > uuqq), "tau channels" ZH — > qqr + r~ and ZH — > r + r~bb (ZH — > r + r~qq), and 
"lepton channels" ZH — > e + e~bb and ZH — > /i + /i~bb (ZH — > e + e~qq and ZH — » ji + n~qq). Those final 
states together cover approximately 80% of the total production rate (see Tabled]) for a Higgs boson 
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Table 4: Fraction of production rates in the different search channels at LEP2. The number outside 
(inside) brackets correspond to H — y qq (H — y bb) decays. 



Channel 


iwH — ou vje v 


/l/f 11C PoV 

iriH — llOVjreV 


TAUT 1 1 nr C 

loui jeis 


u^uu ) /o 


I ^OU J /0 


missing energy 


18(17)% 


16(14)% 


taus 


9(9)% 


6(5)% 


leptons 


6(6)% 


8(8)% 


sum 


97(91)% 


87(77)% 



with a mass of 115 GeV, when only taking into account the decay H — y bb. Given the long time intervals 
between bunch crossings of 22 /is, the instantaneous luminosity of 10 32 cm _2 s _1 , and the size of the total 
cross section, triggering on all signal topologies was not a problem. Due to the presence of irreducible 
backgrounds from pair production of weak gauge bosons, in particular for CM energies of 183 GeV and 
larger, it was not possible to design analyses that yield background expectations at the level of few 
events as at LEP1, while retaining a large fraction of the signal process. Extensive use was made of 
largely improved capabilities to identify b-flavoured jets in order to suppress backgrounds, in particular 
from WW production. The reconstruction of the mass of Higgs-boson candidates was greatly improved 
in many channels by kinematic fits which required energy and momentum conservation and often the 
consistency of a di-fermion mass with M z . Many analyses exploited the advantages of multivariate 
techniques, which were trained for various different Higgs-boson mass hypotheses. In several cases the 
mass- sensitive observables together with the response of the multivariate technique were combined in a 
final discriminant, which was used in the hypothesis tests. 

Below the main characteristics of the four search channels, which are ordered by their sensitivity, 
and the most important means to discriminate them from the dominant backgrounds are discussed: 

• "Four-jet channel": The topology is characterized by four jets, two consistent with the decay Z —y 
qq and two consistent with the decay H — y bb. Requiring two jets to be identified as b-flavoured 
suppresses strongly backgrounds from WW production and partially from ZZ production, as the 
branching ratios for Z — y bb of 15% is much smaller than the corresponding one for the Higgs 
boson. The sensitivity is enhanced by a precise reconstruction of the di-jet invariant mass of 
the jets assigned to the Higgs-boson decay. Exploiting energy and momentum conservation and 
constraining the other di-jet mass to Mz, significantly improves the mass resolution. Contributions 
from e + e~ —y Z/7* —y 4 jets are small for CM energies beyond 183 GeV. 

• "Missing-energy channel" : The topology is characterized by large missing energy, missing (trans- 
verse) momentum, and missing mass (derived from energy and momentum conservation) due to 
the decay Z — y uu, two acoplanar b-flavoured jets from the decay of the Higgs boson, and no 
leptons being present. Requiring b-tagged jets greatly reduces backgrounds from WW —y qqiis, 
where the charged lepton escapes detection, and partially helps to reduce the contribution from 
ZZ —y qqvv. The reconstruction of the mass of the Higgs boson candidate from the di-jet system 
is often significantly improved by exploiting energy and momentum conservation and constraining 
the missing mass to Mz. The remaining background is dominated by di-boson production. 

• "Lepton channels": The topology is characterized by two acoplanar leptons consistent with the 
decay Z — y and two acoplanar b-flavoured jets from the decay of the Higgs boson. Requiring 
b-tagged jets reduces background from ZZ — y qq£ + £~. The mass of the Higgs-boson candidate is 
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Table 5: Excluded mass ranges in the individual experiments at various stages during the LEP2 pro- 
gramme: data set used, observed (expected) mass limit, and reference. The entry "-" indicates that 
this information is not published. 



CM energy 


ALEPH 


DELPHI 


L3 


OPAL 


< 161 GeV 

< 172 GeV 

< 183 GeV 

< 189 GeV 

< 202 GeV 

< 209 GeV 


70.7 (-) 11711 

V / 1 — ■ i 

87.9 (85.3) jl78j 
92.9 (95.9) |182j 
107.7 (107.8) |186j 
111.1 (114.2) [172] 


66.2 (-) 11751 

V / 1 — ■ i 

85.7 (86.5) |l79j 
94.6 (94.4) |183j 
107.3 (106.4) |187j 
114.3 (113.5) [189] 


69.5 (-) 11761 

V / 1 — ■ i 

87.6 (86.8) |180j 
95.3 (94.8) |184j 
107.0 (105.2) |188j 

- [no] 


65.0 (-) [174] 
69.4 (65) 11771 

\ / 1 — ■ i 

88.3 (86.1) [181] 
91.0 (94.9) [185] 

109.7 (112.5 ) [191] 


final 


111.5 (114.2) |192] 


114.1 (113.3) [193] 


112.0 (112.4) [194] 


112.7 (112.7) [195] 



reconstructed as the recoil to the di-lepton system. In some cases kinematic fits exploiting energy 
and momentum resolution improve the sensitivity. For CM energies larger than 183 GeV the final 
background is, to a very large extent, due to ZZ production. 

• "Tau channels" : The topology is characterized by two hadronic jets and two tau-lepton candidates, 
where one pair is assigned to the Z decay and the other to the H decay. The assignment to the 
ZH — > qqr + T~ or ZH — > T + r~bb channels is based on the invariant masses of the tau and jet 
systems, and the quality of kinematic fits requiring energy and momentum conservation and 
applying the Mi constraint. B-tagging improves the sensitivity in the ZH — > T + r~bb final state. 
Selected background events are largely due to di-boson production. 

After all selection requirement a typical signal-to-background ratio of the order between 2-to-l and 
1-to-l could be achieved at the edge of the LEP sensitivity for masses of 115 GeV still retaining 30—60% 
of the eventually produced signal events depending on the final state. No significant deviations from 
the background-only hypothesis were observed in the data collected up to 1999 (see first four rows in 
Table |5]). When comparing the sensitivity of the four experiments one should keep in mind that slightly 
different techniques for deriving the limits are used: for example ALEPH uses the technique discussed in 
Ref. [196j instead of the CL S technique yielding a better expected limit by up to 0.5 GeV. Shortly after 
the end of data taking in 2000, ALEPH |172j reported an excess of events, dominated by the observation 
in the "four-jet channel", with a minimal local p-value of 1.5 x 10~ 3 for Mh = 116 GeV and stated that 
the observed rate would be compatible with Higgs-boson production for Mh = 114 GeV. The other 
three experiments also published their initial findings based on 2000 data with focus on the mass range 
around 115 GeV. L3 |190j reported an excess with a p- value at Mh = 114.5 GeV for the background-only 
hypothesis of 0.09 (for signal+background of 0.62). DELPHI and OPAL did not observe any excess. 
DELPHI [189] preferred the background-only hypothesis (p-value = 0.23) over the signal+background 
hypothesis (p-value = 0.03) for Mh = 115 GeV, and OPAL [191] reported a similar p-value of 0.2 (0.4) 
for the background-only (signal+background) hypothesis. The corresponding observed and expected 
mass limits are shown in the fifth row of Table [5j The final LEP analyses (see last row of Table |5]) used 
the latest and more accurate determinations of the LEP beam energies, better detector calibrations, 
a higher-statistics sample for simulated events and in some cases performed an optimization of the 
selection strategies for high Mh, whereas non-optimized analyses were used in the first analysis of the 
2000 data. 

The final mass distribution obtained from the combination of the four LEP experiments after tight 
selection requirements is shown in Figure [TBI (top-left). The slight excess at large values is dominated 
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Figure 13: Final results of the Higgs-boson searches at LEP [197J: mass distribution of tightly selected 
Higgs-boson candidates compared to signal and background expectations (top- left), mean value of the 
test statistic —2 In Q as function of Mh for data and expectations from the signal+background hypothe- 
sis and the background-only-hypothesis with 68% and 95% CL error bands (top-right), local p- values as 
function of Mh for data and expectations from the background-only hypothesis and signal+background 
hypothesis with 68% and 95% CL error bands (bottom- left), CLs-values as function of Mh for data and 
expectations from the signal+background hypothesis and background-only hypothesis with 68% and 
95% CL error bands (bottom-right). 



by the "four-jet channel" in ALEPH. The expected and observed values of the test statistic —2 In Q 
depending on M H (Figure [T31 top-right), show the limited capabilities to distinguish the two hypothe- 
ses for large Higgs-boson masses. The two bands show the 68% and 95% probability bands for the 
background-only hypothesis. The final local p-values for the combination of all experiments are shown 
in Figure [13] (bottom,left) and detailed in Table [6] for the combination and the individual experiments. 
No significant excess is observed in the combination. The p- values for the two alternative hypotheses are 
very similar, indicating again the difficulty to discriminate the background-only and signal+background 
hypotheses for Mh = 115 GeV. All SM Higgs-boson mass hypotheses below 114.4 GeV were excluded 
with at least 95% CL from the combination of the four experiments (see Figure [T3l bottom-right) [197J. 
The expected combined limit was 115.3 GeV. The very small differences in the limits from the LEP 
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Table 6: Final results of the searches at LEP2: local p-values for the consistency with the background- 
only and signal+background hypotheses assuming Mh = 115 GeV, expected and observed mass limits 
as derived in the LEP Higgs Working Group |197j . 





p 

b-only 


H i i, 

A S-f- D 


exp. limit 


obs. limit 


LEP 


0.09 


0.15 


115 3 GeV 


114.4 GeV 


ALEPH 


3.3 x 10~ 3 


0.87 


113.5 GeV 


111.5 GeV 


DELPHI 


0.79 


0.03 


113.3 GeV 


114.3 GeV 


L3 


0.33 


0.30 


112.4 GeV 


112.0 GeV 


OPAL 


0.50 


0.14 


112.7 GeV 


112.8 GeV 


four jets 


0.05 


0.44 


114.5 GeV 


113.3 GeV 


all but four jets 


0.37 


0.10 


114.2 GeV 


114.2 GeV 



Higgs Working Group for each individual experiment in Table [6] and the final ones published by each 
experiment, as shown in Table El are due to differences in the statistical methods applied. 

At the beginning of the LEP programme no solid limit existed on the mass of the Higgs boson. The 
searches for the SM Higgs boson carried out by the four LEP experiments extended the sensitive range 
well beyond that anticipated at the beginning of the LEP programme. This is due to the higher energy 
achieved and to more sophisticated detectors and analysis techniques. The range below 114.4 GeV was 
and is difficult to probe at past and current hadron colliders. 

6 Higgs-boson production at hadron colliders 

6.1 Higgs-boson production mechanisms and cross-section overview 

The four main production mechanisms for SM Higgs bosons at hadron colliders are illustrated by some 
representative LO diagrams in Figure [TH The size of the respective cross sections depends both on 
the type of colliding hadrons and on the collision energy. Figures [15] and [16] show the total cross 
sections of the various channels for the pp collider Tevatron at its CM energy of \fs = 1.96 TeV and 
for the pp collider LHC at the two energies yfs = 7 TeV and 14 TeV. At the LHC, the energy increase 
from 7 TeV to 8 TeV leads to an increase of 20—30% in the Higgs-boson production cross sections for 
Mh ~ 100— 200 GeV. The transition from 7 TeV to 14 TeV in energy increases the cross sections even 
by a factor of about 3—4 for these Higgs-boson masses, with the exception of ttH production, where 




(a) (b) (c) (d) 

Figure 14: Representative leading-order diagrams for the main SM Higgs-boson production channels 
at hadron colliders, where q and Q denote light and heavy quarks, respectively: (a) gluon fusion, 
(b) Higgs-strahlung, (c) vector-boson fusion, (d) heavy-quark associated production. 
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Figure 15: Cross sections for the various SM Higgs-boson production channels at Tevatron with a CM 
energy of 1.96 TeV, as predicted by the TeV4LHC Section Working Group [198J . 




M H [GeV] M H [GeV] 



Figure 16: Cross sections and respective uncertainties, indicated by the band widths, for the various 
SM Higgs-boson production channels at the LHC with a CM energy of 7 and 14 TeV, as predicted by 
the LHC Higgs Cross Section Working Group |25j. 



the factor is roughly 8. Globally, loop-induced Higgs-boson production via gluon fusion delivers the 
largest cross section owing to the large gluon flux in high-energy proton-(anti)proton collisions. The 
respective cross section is typically an order of magnitude or more larger than the remaining production 
cross sections. At the LHC, the vector-boson fusion (VBF) cross section competes in size with gluon 
fusion only for Higgs bosons as heavy as nearly 1 TeV. VBF delivers the second largest cross section, 
showing a much slower decrease with increasing Higgs-boson mass owing to its t-channel dominance, 
which leads to a logarithmic rise of its partonic cross section with increasing partonic CM energy 



40 



Table 7: Theoretical (THU) and parametric PDF+a s uncertainties (PU) for the total Higgs-boson 
production cross sections at the LHC with CM energies 7TeV and 14TeV, as assessed by the LHC 
Higgs Cross Section Working Group [22] employing the PDF4LHC [21] recipe for PU, as well as the 
typical size of radiative corrections of the strong (QCD) and electroweak (EW) interactions. The colour 
coding NLO/NNLO/NNLO+ refers to the respective perturbative order included in the predictions, 
where NNLO+ means that resummations beyond the fixed-order correction are included. 









LHC @ V 


s = 7TeV 






LHC @ v 


fl = 14 TeV 








uncertainties 


corrections 


uncertainties 


corrections 




M H [GeV] 


THU 


PU 


QCD 


EW 


THU 


PU 


QCD 


EW 


ggF 


< 500 


6-10% 


8-10% 


>100% 


5% 


6-14% 


7% 


>100% 


5% 


VBF 


< 500 


1% 


2-7% 


5% 


5% 


1% 


3-4% 


5% 


5% 


HW 


< 200 


1% 


3-4% 


30% 


5-10% 


1% 


3-4% 


30% 5 


-10% 


HZ 


< 200 


1-2% 


3-4% 


40% 


5% 


2-4% 


3-4% 


45% 


5% 


ttH 


< 200 


10% 


9% 


5% 


? 


10% 


9% 


15-20% 


? 



At Tevatron the Higgs-strahlung channels of HW/HZ production compete with VBF in size for Higgs 
masses Mh J$ 100—200 GeV mainly due to the different combinations of PDFs. For a pp initial state 
high-energy qq collision, which is needed for Higgs-strahlung in LO, is preferred over qq scattering, since 
in this case the qq channel can proceed via two valence quarks, which carry much more momentum 
than sea quarks on average. The total Higgs-strahlung cross section for HW production, where the 
sum over W ± is taken, is larger than the one for HZ production by roughly a factor of two. The 
smallest relevant cross section in all cases is provided by Higgs production in association with tt pairs, 
whose suppression is mainly due to the large invariant mass required to produce the three-particle final 
state of heavy objects. Similar to pure tt pair production, ttH production is largely dominated by qq 
annihilation at Tevatron, but by gg fusion at the LHC. In the transition from Tevatron to the LHC and 
with increasing LHC energy, the suppression of the ttH cross section steadily decreases with respect to 
the other channels, since the phase-space suppression of the heavy ttH final state is fading for larger 
collider energy. The main motivation to measure the ttH cross section clearly rests in the direct access to 
the top Yukawa coupling, without contamination from other couplings; even a qualitative measurement 
at the LHC would be a great success. 

Predictions for hadronic collisions in general involve several serious sources of uncertainties that are 
tied to the hadronic environment, as explained in more detail in Section l2~4l Perturbative corrections, 
especially of the strong interactions, have to be taken into account as much as possible, in order to 
minimize the uncertainties, and the residual uncertainties, which are due to missing higher-order effects 
and parametric errors have to be quantified carefully. Figure [TBI illustrates the error estimate as assessed 
by the LHC Higgs Cross Section Working Group for the total cross sections at the LHC [25], compris- 
ing both theoretical and parametric errors (PDF+a s ). Table [7J gives a brief overview of theoretical 
and parametric PDF+a s uncertainties, for the latter following the recipe of the PDF4LHC Working 
Group [21]; the explicit numbers are based on the results given in Ref. [22]. For later reference, the 
table also illustrates the typical size of the known QCD and EW corrections. A similar, but rather 
conservative estimate of cross-section uncertainties for the Tevatron can be found in Ref. |199j . 

The overview over the total cross sections shown in Figures [15] and [16] can only give a rough idea 
about the importance the respective production channels for Higgs-boson discovery or for later precision 
studies. The total cross sections provide the total production rates for Higgs bosons, but in practice 
the rates of potentially observable events is relevant. Moreover, the signatures left in the detectors are 
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rather different for the various production channels: 

• Single-Higgs production in gg fusion strongly prefers Higgs bosons at low transverse momenta plus 
some hadronic activity recoiling against the Higgs boson. Owing to the overwhelming hadronic 
background it is impossible to access the H — > bb decay channel here. Instead the H — > 77 decay 
yields a promising signal at small Higgs-boson masses in spite of its small branching ratios. At 
higher values of Mh the decays H — > WW/ZZ — > 4fermions via weak gauge bosons lead to rather 
clean signatures, in particular in the "gold-plated" four-muon channel of the ZZ decay. 

• VBF offers the possibility to tag on the two accompanying jets which strongly prefer the forward- 
backward directions owing to the t-channel topology. Colour exchange between the two colliding 
partons is largely suppressed also in higher orders, further pronouncing this signature. Special 
VBF selection cuts on the two jets exploit this feature to suppress background. Apart from the 
WW/ZZ decays, VBF offers the best possibility to search for Higgs-boson decays into t + t~ pairs. 

• Higgs-strahlung represents the main search channel for SM Higgs bosons near the LEP exclusion 
limit at Mh = 114 GeV at Tevatron, but at the LHC the higher background rates pose serious 
problems. In order to make use of this channel, a sophisticated novel search technique is promising 
that first focuses on "fat jets" (with large cones) at large transverse momenta which show an 
internal structure with B mesons [200J. Heavy objects decaying into a bb quark pair, such as a 
light SM Higgs boson, are very likely to fit into this selection pattern. 

• ttH production seems at best to be accessible at very low Higgs-boson masses owing to the cross- 
section suppression, so that one is forced to search for the H — )• bb decay. This leads to the 
partonic final state ttH — > WWbbbb — > 4/ + 4b, where / stands for a four-fermion final state 
of the type 4£, 2£2q, or 4q. It turns out that it will be extremely difficult to establish a signal 
over the serious background created by ttbb and tt + 2 jet production. At present, here the hope 
again rests on a selection focusing on fat jets with a bb substructure |2Ulj . However, owing to 
these difficulties, ttH production cannot contribute to the Higgs-boson discovery, but would be 
very interesting in later precision analyses to assess Higgs-boson couplings. 

• Finally, according to the SM, Higgs-boson production with bb pairs is most likely not accessible 
at hadron colliders owing to the overwhelming hadronic background at hadron colliders. However, 
it plays an important role in extensions of the SM where the Hbb Yukawa coupling is enhanced, 
such as in supersymmetric or non-supersymmetric two-Higgs doublet models. 

6.2 Survey of precision calculations 

The features reviewed above make it clear that theoretical predictions have to account for kinematical 
details of the various production mechanisms and of the Higgs-boson decay channels, because many of 
the specific event selections heavily make use of those properties. 

For the individual Higgs-boson production channels in the SM, theory predictions are accurate 
within 5—20%, depending in detail on the production mechanism: 

• Gluon fusion pp/pp H + X [72ll74ll8lira M[20^226] 

Higgs-boson production via gluon fusion is strongly dominated by top-quark loops in the SM and 
receives very large QCD corrections, which are mainly due to soft-gluon exchange and emission. 
The NLO corrections, which are known with their full quark- mass dependence [T2"|l2"03j. increase 
the cross section by 80—100%. The NNLO QCD corrections, which were first evaluated in the 
large-top- mass limit |204j and later on refined by top-mass corrections |205j . add another 25% to 
the cross section. The still significant scale uncertainty of the NNLO QCD cross section can be 
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reduced by soft-gluon resummation. The resummation effects, known up to the next-to-next-to- 
next-to-leading logarithmic (N 3 LL) level [206] . add another 6—9% to the cross section at the LHC, 
leading to the 6— 10%(6— 14%) of residual scale uncertainties at 7(14) TeV quoted in Table CD The 
seemingly bad convergence in the sequence LO— t-NLO— i-NNLO of fixed-order predictions can be 
cured by an analytical continuation of the cross-section result to an imaginary renormalization 
scale, which results in a resummation of large corrections cx (a s 7r 2 ) n [207j . However, the optimistic 
estimate of uncertainties in the effective- field theory prediction of Ref. [207] . which is also accurate 
to the N 3 LL level, is still under debate [25] . 

EW corrections, which are known at NLO [74 1I8T1 12081 1209] . turn out to be of the size of ~ 5% 
(relative to LO) and strongly depend on the Higgs-boson mass. The extreme size ~ 100% of 
the QCD corrections poses the question how to combine the NLO 0(a s ) QCD and 0(a) EW 
corrections: additively or in factorized form? An investigation of mixed QCD-EW corrections 
for small Mh clearly favours the factorized approach [210] . Electroweak corrections to Higgs+jet 
production [211] influence the total single-Higgs production cross section at the negligible < 1% 
level, but become relevant at large Higgs-boson transverse momenta. 

Finally, at the level of some percent, interference effects between Higgs-boson signal and irreducible 
background, and Higgs-boson off-shell effects should be taken into account, as studied for gg — >• 
H -> WW/ZZ in Refs. [5511511212] . 

Various predictions of total Higgs-boson production cross sections [901 11991 T207] 1210] 1213] were 
described and partially updated in the first review [25] of the LHC Higgs Cross Section Working 
Group, where many details, in particular, about the various error estimates can be found. The 
second review [25] summarizes differential cross sections, such as the Higgs-boson-transverse- 
momentum (pt,h) distribution, and related cut-based observables, such as cross sections with a 
jet veto or definite jet multiplicities — a broad subject where we can only make brief statements. 

Predictions for the j?t,h distribution were first worked out in fixed-order NLO QCD [214] and were 
refined later by soft-gluon resummations [21511216] . For high Pt,h ; the application of the large-top- 
quark-mass limit in predictions becomes subtle, and calculations should include the full quark-mass 
dependence as far as possible; errors from missing corrections due to finite quark masses do not 
necessarily show up in the scale uncertainty. Moreover, differential if-factors [218,219,221], i.e. 
corrections to distributions, in general show significant variations over phase space or depend on 
cuts, clearly disfavouring the use of a simple uniform rescaling of distributions by global .fT-factors. 
Even for integrated quantities, in particular, if vetoes or cuts are involved, a naive scale variation 
can underestimate theoretical uncertainties significantly. For the jet multiplicities and jet vetoes 
the discussion of appropriate error estimates and resummations started in Ref. [26J and went on 
in the literature since then, with proposals for solutions based on traditional resummations [227] 
and effective-theory approaches [217] . 

Several public programmes that include higher-order corrections to Higgs-boson production via 
gluon fusion are used in the experimental anaylses: 

— Higlu [72] for calculating the total NLO QCD cross section with the full mass dependence; 

- MCFM [222J[22I] for Higgs+2jet production at NLO QCD; 

- FEHlP [2TB] and HNNLO J2T9] for differential quantities, based on the NNLO QCD calcu- 
lation in the heavy-top limit, but supporting the full kinematical information on the Higgs- 
boson decays H — > 77 and H — >■ WW/ZZ — » 41eptons; 

- HqT [2T6] and HRes [221] for the NNLL resummation matched to NLO QCD for Higgs- 
transverse-momentum spectrum, where HRes supports the decays H — > 77 and H — > 
WW/ZZ -> 41eptons; 
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— iHixs [55 |I220] for the inclusive Higgs-boson production cross section using state-of-the-art 
fixed-order results, i.e. QCD corrections to NNLO, EW and mixed QCD-EW corrections, as 
well as quark-mass and finite width effects; 

- MCONLO [225] and Powheg [226] for differential predictions in NLO QCD properly 
matched to QCD parton showers. 

• Vector-boson fusion pp/pp — >■ H + 2jets + X [8^ imiT29ll223ll224ll228fM0] 

Higgs-boson production via VBF experimentally means Higgs-boson production in association 
with two hard jets which are forward-backward pointing with some rapidity gap and an appro- 
priate central-jet veto, in order to drastically reduce background [229] . The VBF cuts suppress 
H + 2 jet production via top-quark- loop-induced gg fusion [230] as well, which for M H = 120 GeV 
contributes about 4—5% |23f j to the cross section with a non-negligible residual NLO scale un- 
certainty of about 35% [223 . 1224] . While s-channel contributions from HW/HZ with hadronically 
decaying W/Z bosons deliver up to ~ 30% to the total cross section at low M H ~ 120 GeV, VBF 
cuts suppress this Higgs-strahlung contamination to < 0.6%. Thus, H + 2 jet production with 
VBF cuts is nearly a pure EW process proceeding via W- or Z-boson fusion to a Higgs boson with 
almost no colour exchange between two forward-scattered (anti) quark lines. 

The suppressed colour exchange and the forward kinematics suggest that QCD corrections are 
small after an appropriate scale setting (to small scales of the order of Mw) in the PDF redefinition 
of QCD factorization. In fact this expectation is confirmed by explicit NLO QCD [89,232-234] 
and NNLO QCD [23514237] calculations, which are available for differential cross sections in the 
former case and for the total cross section in the latter. The NLO QCD corrections are of 
the order of 5—10%, depending in detail on the VBF cuts, with interferences between t- and 
w-channels — including gluon exchange between the (anti)quark lines at NLO — at the negligible 
level of < 0.1% p34j. This suppression of colour exchange is also confirmed at NNLO, where large 
parts of the corrections with non-trivial colour flow were calculated or estimated [235 , 1237] . so that 
the NNLO QCD corrections can be treated in the so-called structure-function approach as used in 
Refs. |236tl237j for the total cross section. Likewise, interference effects between amplitudes with 
top-quark-loop-induced ggH couplings and VBF amplitudes with one-loop gluon exchange turn 
out to be negligible as well [238]. The remaining scale uncertainty of the total cross section at 
NNLO QCD, which is at the percent level, turns out to be small as compared to the parametric 
PDF+a s uncertainty of 2— 7%(3— 4%) at the LHC at 7(14) TeV. For differential quantities, i.e. 
for cross sections with VBF cuts or for distributions, QCD corrections at fixed order are only 
known to NLO, but including the matching to QCD parton showers at this order [239J. 

The NLO EW corrections have been calculated for all t/u- and s-channel contributions, including 
all interferences, in Ref. |234] and for the t/u-channels in Ref. |240] . where both calculations are 
valid for integrated and differential quantities. In the presence of VBF cuts, the approximation 
by t/w-channels only is sufficient. The EW corrections to integrated cross sections amount to 
5—10% and are, thus, of the same generic size as the QCD corrections, but in distributions at 
high transverse momenta of the Higgs boson or the jets the EW corrections can even be larger. 
Similar to the situation in the related decays H — > WW/ZZ — > 4f, discussed in Section ¥Z15\ the 
leading two-loop corrections of the order G^M H [ST] become of the same size as their one-loop 
counterpart of order G M M H at Mh ~ 700 GeV, signalling the breakdown of perturbation theory. 
For Mh ^ 500 GeV, on the other hand, perturbation theory is safely applicable. 

The described state-of-the-art is widely available by public programmes. The fixed-order NLO 
QCD+EW predictions can be obtained with the Monte Carlo programmes VBF@NLO |233[ 
1240(1241] and HAWK |234[l242j . the former in the t/w-channel approximation, the latter including 
s-channels and all interferences as well. The (approximated) NNLO QCD corrections to the 
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total cross section can be obtained with the programme VBF@NNLO [236, 1237(1243] . Presently, 
the best prediction for differential distributions should be obtained upon reweighting the parton- 
shower-improved NLO predictions of Powheg |239] with the relative EW correction from Hawk 
or VBFONLO. 

• Higgs-strahlung pp/pp HW/HZ + X [213H253] 

As far as higher-order QCD corrections are concerned, Higgs-strahlung is very similar to the 
well-known Drell-Yan process. In LO there is just quark-antiquark annihilation to an off-shell 
weak gauge boson, qq ->■ W*/Z* ->■ HW/HZ, and the structure and size (~ 30%) of NLO QCD 
corrections [245 j are identical to the Drell-Yan case. At NNLO QCD, the Drell-Yan-like correc- 
tions, which are known for the total HW/HZ cross sections [246] and the differential HW cross 
section |247] , dominate, but there are additional corrections to Higgs-strahlung that have no coun- 
terparts in Drell-Yan production. Firstly, loop-induced gg fusion contributes to HZ production 
about 2-696(4-12%) for CM energies of 7TeV(14TeV) at the LHC |23Sfi and secondly dia- 
grams involving Higgs couplings to top-quark loops correct the HW/HZ cross sections by another 
1—3% |249] . As far as EW corrections are concerned, Higgs-strahlung is more involved than the 
Drell-Yan process. The NLO EW corrections to the total cross sections have been calculated in 
Ref. [250J and are of the order of —(5—10)%, only weakly depending on the collider type and 
energy. The NLO EW corrections to differential cross sections, calculated in Ref. [251], can be 
significantly larger, in particular for the kinematics relevant in the "boosted-Higgs analysis" which 
requires large transverse Higgs momenta. For instance, the EW corrections to HW production 
with Mn = 120 GeV are 14% for p T ,H > 200 GeV. 

From Table [7] one can read off that theoretical uncertainties at the LHC are of the order of 
l-2%(2-4%) and the PDF+a s uncertainties of the order of 3-4% for a CM energy of 7(14) TeV 
for the total HW/HZ cross sections (see Ref. [25] for details). Note, however, that the uncertainties 
for the scenario of a boosted Higgs boson (pr,H ~ 200 GeV) are somewhat larger [26J. Specifically, 
the issue of final-state radiation off the Higgs-boson decay products deserves particular attention 
in the boosted-Higgs scenario [252J. The PDF uncertainties are estimated to ~ 5%, and for HZ 
production, where differential QCD corrections are only known to NLO, the theoretical uncertainty 
is estimated to ~ 7%; a reduction of these theoretical uncertainties to the few-percent level requires 
the inclusion of NNLO QCD corrections. 

To a large extent, the described precision calculations are available via public programmes. The 
NLO QCD corrections for total and differential cross sections can be calculated with V2HV [254] 
and MCFM [222], respectively, and the total cross section, including NNLO QCD and NLO EW 
corrections, with the programme VH@NNLO [253] . Finally, the NLO QCD+EW corrections to 
fully differential observables (with decaying W/Z bosons) can be computed with the Monte Carlo 
programme HAWK [242] . 

• Associated Higgs-boson production with tt pairs pp/pp — > ttH + X [201U255H 258J 

Both the total and differential cross sections to ttH production are known to NLO QCD [256J. 
The size of the QCD corrections is moderate, i.e. typically 10—20% (positive at the LHC and 
negative at Tevatron), which holds also true for differential distributions if appropriate dynamical 
scales are used. Theoretical uncertainties from missing higher-order effects as well as PDF+a; s 
uncertainties are estimated to be of the order of ~ 10% for total cross sections at the LHC (see 
Table [7] and Ref. [25J). For differential distributions the theoretical uncertainty is somewhat larger 
(10—20%) and together with the parametric PDF+a s uncertainties can add up to 20—50% [26] . 

8 Thcsc numbers, which were used also in Figure [THJ are based on the LO one-loop diagrams of the gg channel. A recent 
evaluation [248] of the two-loop NLO corrections in the heavy-top limit, however, roughly doubles this contribution, i.e. 
in order to reach the precision tag of Table [7] these corrections should be taken into account. 
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The NLO QCD predictions were matched to parton showers, employing the MC@NLO |257j and 
Powheg [258] concepts. The results obtained with the two shower variants typically agree to 
better than 10%. Electroweak corrections to ttH production are yet unknown, but their effect 
should be covered by the size of the quoted QCD uncertainties. 

As already mentioned above, suppressing the background to ttH production is a real challenge, 
since the final state involves four b quarks, where three or four bottom quarks have to be tagged, 
and two out of the four have to be identified as H — > bb decay products. To properly simulate 
this theoretically, precise predictions are required that take into account the decays of both the 
top quarks and the Higgs boson. At present, such predictions exist only at LO. Similar to the 
HW/HZ analysis, searching first for fat jets at high transverse momentum and a bb substructure 
can lift the signal-to-background ratio to S/B ~ 1/2.4, assuming three b tags, Mr = 120 GeV, 
and an integrated luminosity of lOOfb -1 , as shown in the parton- level study of Ref. |201j . The 
main background results from ttbb and tt+2jet production within pure QCD, whose cross sections 
have recently been calculated at NLO in Refs. [259] and |260] , respectively. However, the residual 
scale dependence of the background cross-section normalization is about ~ 20%, and a careful 
study of the uncertainties of possible extrapolation procedure from signal-free control regions to 
signal regions has not yet been carried out. To this end, for the background the differential 
information on the top-quark decays has to be taken into account as well. The feasibility of the 
ttH measurement is still not established. 

6.3 Towards couplings analyses after discovery 

Although the difficult hadronic environment of the LHC limits the precision of Higgs-boson studies, 
the global pattern of Higgs couplings could be tested at a qualitative level in a global fit to various 
channels with only mild theory assumptions. Assuming, for instance, that the Higgs-boson decay 
widths into weak gauge bosons are at most as large as their SM values, r H _j. ww / zz < r^ ww y zz , 
which is the case in many models with extended Higgs sectors, the Higgs couplings of a not too heavy 
Higgs boson (110 GeV < M H < 190 GeV) to top quarks, r-leptons, and W and Z bosons could be 
determined to 10—40% in an LHC run at -/s = 14TeV with an integrated luminosity of 300 ftT 1 [ 26 1L 
263J. An increase of the integrated luminosity to 3000 fb _1 reduces the experimental uncertainties in 
the signal strengths a /osm by roughly a factor of 2 [262J, so that the uncertainties of the Higgs-boson 
couplings should be reduced by a factor of about 1.5, since the couplings enter the signal strengths 
quadratically. Constraints on Higgs couplings might be even possible with lower luminosity or tighter 
constraints with the same luminosity if more and more sophisticated analysis techniques are exploited 
(see, e.g., Ref. [2M] ). Currently, the ATLAS [[265] and CMS [266] collaborations follow the strategy of 
a first model-independent analysis of Higgs couplings, which has been put forward by the LHC Higgs 
Cross Section Working Group [267] and is based on a simple rescaling of SM cross sections. Ultimate 
precision Higgs coupling analyses, however, should go beyond such simple phenomenological studies, 
since their informative value is rather limited in the interesting case of significant discrepancies between 
SM and data. Proper model-independent non-standard coupling analyses should be thoroughly based 
on effective-field theories, including higher-order effects in a consistent way as much as possible. First 
steps into this direction are proposed in Ref. [268J. 

Among the various Higgs couplings, the self-interactions are of particular interest, since they deter- 
mine the Higgs potential that drives EW symmetry breaking. As already explained for e + e~ colliders, 
such studies require multiple-Higgs-boson production, e.g. Higgs-boson pair production to access the 
triple-Higgs coupling. The various channels for Higgs-boson pairs at hadron colliders are illustrated 
by some LO diagrams in Figure [T71 The cross sections for (a) gluon fusion [269(1270] . (b) Higgs- 
strahlung [271( 1272]. (c) vector-boson fusion [273(1274] are, however, very small and further suppressed 
by branching ratios of the many-particle final states. At present, it seems that only gluon fusion may 
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(a) (b) (c) 

Figure 17: Representative LO diagrams for the main SM Higgs-boson pair production channels at 
hadron colliders, where q and Q denote light and heavy quarks, respectively: (a) gluon fusion, (b) Higgs- 
strahlung, (c) vector-boson fusion. 

be able to establish a signal for Higgs-boson pairs and to allow for a qualitative measurement of the 
triple-Higgs coupling. Higher-order QCD corrections to gg — > HH are very similar to single-Higgs pro- 
duction gg — >• H, though kinematically more complicated; the NLO QCD corrections to gluon fusion 
to Higgs-boson pairs, which are known in the heavy-top limit |270j . confirm this expectation and yield 
a i^-factor of the order of 2. Based on the HWW decays, a non-vanishing H 3 coupling could be es- 
tablished for Mh = 150—200 GeV after collecting a high integrated luminosity of 300fb _1 at the LHC 
with design the energy of 14TeV [275J. Although the theoretical simulation |27b'j suggests that for a 
low mass of M H ~ 120 GeV the rare decay HH — > bb77 could establish the coupling for an integrated 
luminosity of 600fb~ 1 , an experimental sensitivity study [262] for an integrated luminosity of 3000fb _1 
at a CM energy of 14 TeV indicates that more than one channel will be necessary to empirically establish 
a non- vanishing triple-Higgs-boson coupling. 

Since triple-Higgs-boson production cross sections are way too small |277j at hadron colliders as 
well, an experimental test of quartic Higgs couplings at colliders is out of reach. 



7 Searches at Tevatron 

The Tevatron Run II started in 2001 with a pilot run at a new CM energy of 1.96 TeV. From spring 
2002 until its shutdown in autumn 2011 it provided proton-antiproton collisions corresponding to a data 
set of up to 10 fb" 1 to the two experiments CDF and DO. In contrast to the searches at LEP, where 
the a priori signal-to-background ratio was at the level of 10 -4 to 10 -6 as both signal and background 
processes are mediated by the electroweak interaction, the huge total inelastic cross section due to 
processes mediated by the strong force is more than ten orders of magnitude higher than the cross 
section for Higgs-boson production (see Figure |6]). Owing to the large cross section, the instantaneous 
luminosity of up to 4 x 10 32 cm' 2 s _1 and the short time interval between bunch crossings of 396 ns not 
all interactions can be recorded. In order to limit the recorded data volume to a manageable size severe 
requirements on the presence and quality of final-state objects have to be imposed already at the trigger 
level. This online selection for Higgs-boson searches relies mostly on the presence of a colour-neutral 
object in the final state, i.e. photons, electrons, muons, or significant missing transverse momentum 
arising from neutrinos. Hence, although Higgs-boson topologies with fully hadronic final states, which 
comprise production in gluon fusion with subsequent decays via H — >■ bb or H — >• WW(ZZ) — > 4 quarks, 
lead to the highest production rates, these cannot be exploited at Tevatron. Exceptions to this rule of 
thumb are the searches in the associated production WH(ZH) — > qqbb and in vector-boson fusion with 
the decay H — > bb performed by CDF |278j . where events were recorded by requiring multi-jets and a 
large scalar sum of transverse energies at the trigger level. 
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M H (GeV) 

Figure 18: Production cross sections in the most sensitive search channels at the Tevatron. The numbers 
are extracted from Ref. [279J. 

7.1 Overview of search channels 

Generically, the Higgs-boson mass range 100—200 GeV was investigated by both experiments CDF 
and DO. An overview of the production rates for most of the search channels considered at Tevatron 
is shown in Figure [18j For Mh < 130 GeV the most abundant and most sensitive search channels 
are H — > bb produced in association with a weak gauge boson W or Z, which decay leptonically via 
W — > £v, Z — > £ + £~ , or Z — > vv. Here and in the following £ denotes an electron or muon. The relative 
magnitude of the rates is due to different couplings of the Higgs bosons to W and Z and the branching 
ratios BR(W Iv) = 22%, BR(Z £+£') = 6.7%, and BR(Z -* vv) = 20%. For M H in the range 
130—200 GeV the sensitivity is dominated by the search for H — > W + W~ — > £ + v£~v produced in gluon 
fusion. The discussion below focuses on these four final states. In addition to those high sensitivity 
channels a plenitude of other channels was analyzed which yield cross-section limits for Higgs-boson 
production with a value of at best 7 times the SM prediction, and a sensitivity of a single channel 
(expressed in expected limits) of at least a factor 5 worse compared to a combination of the four channels 
mentioned above. Those channels include H ->- 77 (CDF [280] and DO [281]), H TIL ->• £+£' £'+£'- 
(CDF [282]), qqbb from WH(ZH) and in VBF with the decay H ^ bb (CDF [278]), ttH with H 
bb (CDF [383]), H T+T- (CDF [281], DO [285]), H ->• W+W" ->• £vqq (DO [286]) and W(Z)H -> 
W(Z)WW -> £ ± £ ± + X (CDF [2S7], and DO [255] ). Additional channels are included in the preliminary 
combinations in each experiment and the Tevatron combination. All those channels contribute to 
the overall sensitivity at the level of 10—20%. In most analyses the final discriminant is given by 
an event classifier based on the combination of multivariate techniques, which combines information 
from kinematical properties, such as observables sensitive to the Higgs-boson mass and topological 
information from flavour tagging. Usually the pre-selected event samples, where the presence of the basic 
physics objects of the final state under consideration are required, are divided in categories based on jet 
multiplicity, lepton flavour, multiplicity of b-tagged jets and the tightness of their selection requirements, 
and on additional kinematical properties, such as low or large missing transverse momentum. Often 
the selection strategy using multivariate techniques is a multi-step procedure. In the first step either 
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dominant reducible backgrounds are greatly reduced, or the events are divided in categories, in which a 
specific background is enriched or depleted. In the next selection step the signal can then be separated 
from the remaining backgrounds in a more efficient way. Sometimes a sequence of multivariate event 
classifiers is also applied, which one after the other try to suppress the dominant background classes. The 
final discriminants are then constructed from the responses of various multivariate techniques. This 
sophisticated selection procedure improves significantly the signal-to-background ratios in particular 
categories with respect to a more inclusive analysis and hence the overall sensitivity of the combined 
search. Higgs-boson production is simulated with Pythia [9TJI92], to which a reweighting of the 
transverse-momentum spectrum of the Higgs boson is applied. QCD multi-jet production is estimated 
entirely using daten-driven techniques. Other backgrounds are estimated from simulated events, which 
are generated with a variety of event generators such as Alpgen |289j . Comphep [290] . Herwig |291[ 
[292] . MCONLO [225l[293], Pythia [HUES], and Sherpa [291]. The misidentification probabilities for 
leptons and b-flavoured jets, as measured in data, are applied to the simulation. 

7.2 Searches for H — > bb 

The searches for H — > bb in the mass range 100—150 GeV are performed in three non-overlapping final 
state topologies: (a) 1 lepton and large missing transverse momentum designed for WH — > £vbb, but 
also selecting contributions from ZH — > £ + £~bb, where one lepton is not reconstructed and identified, 

(b) 2 leptons with an invariant mass consistent with Mi optimized for ZH — > £ + £~bb production, and 

(c) leptons and large missing transverse momentum designed for ZH — > vvbb, but also receiving 
contributions from WH — > £vbb and ZH — > £ + £~bb, where one or two leptons are not reconstructed and 
identified. The suppression of background processes, in particular from W(Z) + light-flavour production, 
depends crucially on the performance of the algorithm to identify the b-flavoured jets from the Higgs- 
boson decay. Significant improvements were achieved over the last years, which greatly enhanced the 
sensitivity. The mass of the Higgs-boson candidate can be reconstructed from the invariant mass of the 
two jets assigned to the Higgs-boson decay with a mass resolution of 10—15%. This di-jet mass is an 
important observable to discriminate the signal process from background processes, in particular from 
di-boson production VZ with Z — > bb. In addition to the multiplicity of leptons and the presence of 
missing transverse momentum exactly two or three jets are required with at least one being identified as 
b-flavoured. The events are further divided into categories depending on lepton flavour, lepton quality, 
b-tag multiplicity, and b-tag quality, in order to improve the signal-to-background ratio. For each 
category and Higgs-boson-mass hypothesis multivariate event classifiers are trained, which use the di-jet 
invariant mass and other kinematical quantities as input observables. The response of these multivariate 
techniques is used as the final discriminant. For the purpose of visualization of the sensitivity and results, 
especially when several channels and categories are combined, the bins of each final discriminating 
observable are ordered with respect to the expected signal-to-background (s/b) ratio. The expected 
and observed event yields in each ln(s/b) bin are then added over the discriminating variables in all 
non-overlapping final states and event categories. 

The expected and observed limits on the signal-strength parameter /i = <t/(t sm (called R in the 
Tevatron analyses), which relates the excluded cross section o to the SM expectation o"sm for Higgs- 
boson production, from the analysis of full data set collected during Run II at Tevatron (9.45 fb _1 in 
CDF and 9.5— 9.7 fb -1 in DO) are summarized in Table [8] In both experiments and in all three search 
topologies the observed limit is weaker than the expected one over a broad mass range. The deviation 
is less than 1 standard deviations over the full mass range in the DO £vbb |296] and vvbb searches [298J , 
and within 2 standard deviations in the DO £ + £~bb |300] search and the CDF £vbb [295J and vvbb [299J 
searches. In the CDF £ + £~bb |299] search a maximum excess corresponding to a significance of 2.4 is 
observed at Mh= 135 GeV, reduced to 2.1 when taking into account the "look-elsewhere effect" (LEE). 
When combining the three channels, CDF [302J observes a limit weaker than the expected by 2 standard 
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Table 8: Expected and observed limits on the signal-strength parameter fx = u/usm m the search for 
H — > bb in DO and CDF. The entry "— " indicates that this information is not available. 



channel 


CDF 
obs. (exp.) limit 

at Af TT — 11 K /I Ot; HpV 
ell iV-Z H — XltJ/l^tJ vjtg V 


DO 

obs. (exp.) limit 

at -^-*H — -L / -L^<J VJC V 


WH (H bb, Z -»■ Iv) 
ZH (H ->■ bb, Z i/i/) 
ZH (H ->■ bb, Z -> 


3.1 (2.0) / 4.9 (2.8) [295] 
2.7 (2.7) / 6.7 (3.6) |297j 
4.7 (2.7) / 7.1 (3.9) [2SS] 


3.7 (3.2) / 5.2 (4.7) |296j 
3.0 (2.7) / 4.3 (3.9) |298j 
4.3 (3.7) / 7.1 (5.1) [300] 


combined 


-/- 


1.9 (1.6) / 3.2 (2.3) [301] 



deviations for Mh = 110— 150 GeV and a maximal deviation from the background-only hypothesis for 
Mh = 135 GeV with a significance of 2.7 (2.5 after including the LEE in the mass range 115—150). 
In DO [301J the combination yields a limit weaker than the expected by 1—1.7 standard deviations 
in the mass range 120—145 GeV and a maximal deviation from the background-only hypothesis for 
Mh = 135 GeV with a significance of 1.7 (1.5 with LEE in the mass range 115—150 GeV). DO excludes 
the mass range Mh = 100— 102 GeV. Both experiments conclude that the excess is compatible with 
the production of the SM Higgs boson with a mass of 125 GeV within large uncertainties. The best fit 
values for the signal strength at Mh = 125 GeV is fx = l-2l}j in DO, and in CDF fx is approximately 
2.5 consistent with the SM expectation of 1 within 2 standard deviations. 

The combination of the observed event yields from CDF and DO |303] in the ln(s/b) distribution of 
H — > bb searches after subtraction of the expected background in Figure [TH] (top-left) shows a positive 
deviation from the expectation of zero for the background-only hypothesis in the bins of large signal- 
to-background ratio. The observed limit on the signal strength fx is weaker than the expected by more 
than two standard deviation in the range 115— 145 GeV, as shown in Figure [19] (top-right). Under 
the background-only hypothesis one expects to exclude the Higgs-boson-mass hypothesis in the range 
100—116 GeV. Assuming a Higgs boson with mass of 125 GeV no excluded range is expected. The data 
allow the mass hypothesis below 106 GeV to be excluded. The local p- values in testing the background- 
only hypothesis as a function of Mh (see Figure [T9l top- left) show deviations from the background-only 
hypothesis with a significance of larger than 2 in the mass range 115— 145 GeV with the smallest p- 
value at 135 GeV corresponding to a local significance of 3.3 (2.8 at Mh = 125 GeV). After taking 
into account the LEE in the mass range 115— 150 GeV the significance is decreased to 3.1. The best 
fit value for the signal cross section obtained in data is compared to various types of expectations for 
signal+background hypothesis in Figure [T9l (bottom-right) . The fit values agree with the expectation for 
the SM Higgs boson at the 68 (95)% CL in the approximate mass ranges 112— 120 GeV (105— 132 GeV 
and 140— 150 GeV). As M H = 135 GeV was already excluded by both LHC experiments and guided 
by the excess observed in each experiment at the LHC with a local significance of approximately 3 
at Mh ~ 125 GeV already in the analysis of the 2011 data, the Tevatron experiments quote a best 
fit value for the signal cross section (cx W h + ^zh) x BR(H bb) at M H = 125 GeV of 0.23lg;{jgpb 
to be compared to the SM Higgs boson of 0.12 ± 0.01 pb, consistent within 1.5 standard deviations. 
Despite having no sensitivity to exclude the SM Higgs boson hypothesis for masses larger than 116 GeV, 
the Tevatron observes a broad excess in the H — > bb final state, consistent with a wide range of SM 
Higgs-boson-mass hypothesis. Based on these findings the Tevatron experiments reported evidence for 
the production of a new particle produced in association with weak bosons and decaying to a b-quark 
pair [303J. 
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Figure 19: Results of the search for H — > bb at Tevatron from the combination of CDF and DO full 
Run II data [303] . Expected and observed event yields classified by signal-to-background ratio after 
background subtraction (top-left), exclusion limits on the signal strength R = fi = a/a^m as function 
of Mh (top- right), local p- value for the background-only hypothesis as function of Mh (bottom- left), 
and best fit values for the signal cross section as a function of Mh, compared to the expectation for a 
SM Higgs boson as a function of M H and the expectations for a SM Higgs boson with mass of 125 GeV 
evaluated at different M H (bottom-right). 



7.3 Searches for H W+W~ l+vt~v 

The final-state topology is characterized by two isolated leptons of opposite charge and large missing 
transverse momentum. Due to the presence of two neutrinos in the final state the mass of the Higgs- 
boson candidates cannot be reconstructed and only the transverse mass of the di-lepton and missing- 
transverse-momentum system can be reconstructed. Observables that allow to suppress background 
processes comprise, among others, the di-lepton invariant mass, the transverse mass, and the azimuthal 
opening angle between the two leptons |304j . The latter observables exploit the scalar nature of the 
Higgs boson and the parity violating V — A structure of the W — > Iv decay. Due to the spin correlations 
the signal accumulates at small values of the opening angle and di-lepton invariant mass. The search 
is optimized in different categories, distinguished, e.g., by lepton flavour and jet multiplicity (0, 1, 
>2). The latter aims to enhance the sensitivity for associated WH(ZH) and VBF production with two 
quarks in the final state. The composition of the remaining background is very different in the different 
jet categories: In the (2) jet category the dominating background is given by WW (tt) production. 
The final discriminating observables based on multivariate techniques in each event topology, which are 
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Figure 20: Cross-section limits with respect to the SM expectation for H — > W + W — > £ + v£ v produc- 
tion in CDF |287] (left), in DO |306j (middle), and in the first Tevatron combination [305J (right). 



separately optimized for all Higgs-boson-mass hypotheses tested, are used as input for the statistical 
interpretation of the results. In addition searches for final states with two leptons with same electric 
charge or with three leptons were performed to enhance the sensitivity to WH(ZH) production, where 
the associated W(Z) boson decays to one (two) charged leptons. Due to the comparably large production 
rate, H — > W + W~ — > £ + v£~i? was the first channel to reach the sensitivity to exclude SM Higgs-boson 
cross sections in the vicinity of Mh = 160 GeV. The only published combined Tevatron result [305J 
including H — > W + W~ — > £ + v£~v search channels from spring 2010, which is based on a data sets of 
up to 4.8fb~ 1 in CDF [287] and 5.4fb _1 in DO [306], excludes M H = 162-166 GeV with an expected 
exclusion of Mh — 159—169 GeV (see Figure |20| right). This was the first time that Tevatron expected 
to exclude a Higgs-boson mass range. The individual experiments were not yet able to exclude a mass 
range, but gave observed (expected) limits on the signal-strength parameter at Mh = 165 GeV of 1.29 
(1.20) in CDF (see Figure l20| left) and 1.55 (1.36) in DO (see Figure 1201 middle), to be compared 
with the Tevatron combination 0.93 (0.87). Since then, many preliminary updates of searches for 
H — > W + W~ — > £ + v£rv were performed extending the sensitivity for SM Higgs-boson production to 
a mass window of size of roughly 20 GeV around 2Mw- Only DO published a new result based on a 
data set of 8.6 fb -1 |307] yielding an expected excluded mass range of M H = 159— 169 GeV. However, 
no mass range could be excluded due to an upwards fluctuation of the observed event yields by one 
standard deviation with respect to the background expectation. 



7.4 Preliminary combined results 

Over the years the findings in the plenitude of the decay channels and final-state topologies were 
combined regularly by CDF and DO separately, and also in a common effort as the Tevatron combi- 
nation. The correlations in the systematic uncertainties across channels and also across experiments 
were carefully taken into account. In the following the latest preliminary Tevatron combination (July 
2012) |279] for the mass range Mh = 100—200 GeV is briefly discussed. A comparison of the ex- 
pected and observed event yields classified by signal-to-background ratio after background subtraction 
for Mh = 125 GeV and Mh = 165 GeV is shown in Figure 1211 (top row). The Tevatron combination 
extends the exclusion sensitivity significantly with respect to the one in the individual experiments. 
The mass regions M H = 100— 120 GeV and M H = 139— 184 GeV are expected to be excluded as shown 
in Figure [2T1 (middle- left) . The observed excluded mass ranges are significantly weaker and extend only 
over the mass ranges Mh = 100— 103 GeV and Mh = 147— 180 GeV. The largest deviation from the 
background-only hypothesis is observed for Mh = 120 GeV with a local p-value of 1.5 x 10~ 3 , corre- 
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Figure 21: Preliminary results from the Tevatron combination in July 2012 |279j . Expected and observed 
event yields classified by the signal-to-background ratio after background subtraction for M H = 125 GeV 
(top-left) and M H = 165 GeV (top-right); exclusion limit on the signal strength ct/ctsm as a function of 
Mh (middle-left), local p- value for the background-only hypothesis as a function of Mh (middle-right), 
best fit values of the signal strength as a function Mh (bottom- left), and best fit values for the signal 
strength in different decay modes for Mh = 125 GeV (bottom-right). 
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sponding to a local significance of 3.0 as shown in Figure |2T1 (middle-right) . The excess is dominated by 
the search for H — > bb discussed above. Due to the limited mass resolution the significance is larger than 
2 in the approximate mass range Mh — 115— 145GeV. Hence the influence of the LEE is also rather 
modest with a trial factor of approximately 2 yielding a global significance of 2.5. Another deviation 
from the background-only hypothesis with a significance of larger than 2 is observed for a Higgs-boson 
mass larger than 195 GeV, which was already excluded by the LHC experiments from 2011 data. The 
best fit value for the signal strength as a function of Mh, shown in Figure [211 (bottom- left) . is consistent 
at 68% CL with the expectation for the SM Higgs boson in the mass ranges Mh = 110— 140 GeV and 
approximately Mh = 190—200 GeV. Finally, the best signal-strength parameters in the different decay 
modes H — > bb, H — > WW, and H — > 77 were determined for various fixed values of Mh. The result for 
M H = 125 GeV is shown in Figure |2~T1 (bottom- left). 

8 Searches at the LHC 

The design parameters of the LHC and the two multipurpose experiments ATLAS and CMS are chosen 
in such a way that the whole unexplored mass range, from the limit obtained at LEP up to the upper 
border of < 1 TeV demanded by unitarity (see Section [275]) . can be covered. On 30 March 2010, the 
LHC took first proton-proton collision data at the CM energy of 7 TeV, an energy unprecedented 
at accelerators at that time. Fully hadronic final states, arising from production in gluon fusion or 
weak- vector-boson fusion with subsequent decay via H — )• bb or H — > WW(ZZ) — > 4 quarks yielding 
the largest production rates, are not considered in the search. The overwhelming background and the 
limited trigger capabilities for these final states do not allow for performing a sensitive search. The 
large instantaneous luminosities of up to 6.8 x 10 33 cm _2 s~ 1 (3.7 x 10 33 cm" 2 s~ 1 ) reached in 2012 up 
to June (2011) together with the large total inelastic cross section and the interaction rate of 20 MHz 
lead to an average of 20 (10) overlayed proton-proton interactions in 2012 up to June (2011) often 
referred to as "pile-up". This challenge, which was not of relevance at Tevatron, imposes additional 
challenges on the event reconstruction especially on isolation criteria and a precise determination of the 
missing transverse momentum. Only final states that contain at least one photon, electron, muon, or 
a hadronic tau-lepton decay in association with large missing transverse momentum are considered in 
the searches. The production rates for these topologies at the LHC and the ratio of cross sections at 
different CM energies are shown in Figure [22j In contrast to the searches at Tevatron channels with low 
branching ratios, such as H — > 77 and H — > TIL — > £ + £~£' + £'~ with branching ratios of 0.2% and 0.013% 
for a Higgs-boson mass of 125 GeV are also accessible due to the larger production cross section at the 
LHC. In fact these two channels provide the highest sensitivity over a large mass range. 

The signal cross sections, the branching ratios, and their uncertainties are computed, as detailed in 
Sections 12.61 and |6| following the recommendations of the LHC Higgs Cross Section Working Group [251 
[26] . Simulated events are generated using Powheg interfaced to Pythia for the dominant production 
mechanisms in gluon fusion [226] and in weak- vector-boson fusion [239] . The associated production with 
a weak gauge boson and a pair of top quarks is mainly simulated with Pythia (CMS uses Powheg 
for the decay channel H — > bb). 

The dominant background processes, in particular the irreducible ones apart from a very few ex- 
ceptions, are estimated at least partially with data-driven techniques. Only in some analyses the 
normalization and shape of the final discriminating distribution are completely estimated from data. 
Quite often the normalization is estimated in a control region, defined, e.g., by inverting some of the 
selection criteria, and the prediction for the ratio of the background expectations in the signal and 
control regions as well as the shape of the final discriminant are estimated from simulated events. The 
most advanced event generators are used to get as precise as possible predictions. Those include Alp- 
gen [289], Madgraph [309], MCONLO [225j[293], and Sherpa [291], which perform a matching of 
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M H [GeV] M H [GeV] 

Figure 22: Higgs-boson production rates at the LHC for CM energies of 7TeV in final states which 
are used to search for the Higgs boson (left) [25] and ratio of cross sections at different CM energies 
(right) [508] . 

matrix-element calculations and parton showers at NLO or LO in QCD, and dedicated programmes for 
gluon-induced di-boson production, such as GG2WW [310] and GG2ZZ [311] at LO. For cross checks 
and simulation of the parton shower, of the underlying event, and of pile-up the multipurpose event 
generators Pythia [SUES] and HERWIG [29TIE92] are used. 

8.1 Search with data collected in 2010 

The highest cross sections with final states involving at least one colour-neutral object are obtained 
in the mass range between 150 GeV and 170 GeV in the decay chains H — > W + W~ — > Ivqq and 
H — > W + W~ — > £ + v£~v (see l.h.s. of Figure 122]) . The semi-leptonic decay of the pair of W bosons 
is not usable in this low-mass interval due to the large background from W+jets production. The first 
sensitivity at the LHC hence stemmed from the H — > W + W~ — > £ + v£~v final state. Both experiments 
performed a search in this final state with the data collected in 2010 at a CM energy of 7 TeV correspond- 
ing to integrated luminosities of 36 pb -1 [312J and 35pb _1 |313j . The expected sensitivity in the mass 
range 150—170 GeV was at the level of 2—3 times the SM prediction (see Figj23| middle and right). No 
significant deviation from the expected event yields from background processes was observed. ATLAS 
already investigated several other final states based on 2010 data: H — > 77, H — > ZZ — > £ + £~£ l+ £'~ , 
H — > ZZ — > £ + £~vv, H — 7- ZZ — > £ + £~qq, and H — > W + W~ — > £vqq. The individual sensitivities of each 
of these searches were at the best on the level of 15 times the SM prediction (Fig. I22J left). 

8.2 Searches with data collected in 2011 and up to June in 2012 

In 2011 each of the two experiments collected a data set corresponding to integrated luminosities 
of up to 5.1 fb" 1 at a CM energy of 7 TeV. In 2012 until June a data set corresponding to up to 
5.9 fb -1 per experiment was recorded at a CM energy of 8 TeV. An overview of the search channels 
considered, the considered mass range in each channel, and the integrated luminosities analyzed is 
given in Table |9] Three mass regions can be roughly distinguished: (a) 110 GeV < M H < 150 GeV, (b) 
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Figure 23: Observed and expected exclusion limits on the signal-strength parameter \i obtained from 
2010 collision data: individual channels (left), their combination (middle) in ATLAS [312j . and the 
combination in CMS (right) |313j . 

Table 9: Search channels, mass ranges, integrated luminosities considered by ATLAS and CMS, and 
corresponding references. The results used in the combination of the 2011 data are documented in 
Ref. [314j for ATLAS and in Ref. [315j for CMS. An entry "— " indicates that this search channel was 
not considered for this data set. In 2012 CMS restricted the mass range for all published searches to 
values below 160 GeV. When two mass ranges are given, the search was performed in two disconnected 
ranges. 





ATLAS 


CMS 


channel 


mass range 


/Ldtffb- 1 ] 


mass range 


/Ldt[fb" 1 ] 




M H [GeV] 


2011 


2012 


M H [GeV] 


2011 


2012 


H — > 77 


110-150 


4.8 (316) 


5.9 |317j 


110-150 


5.1 [318] 


5.3 (319] 


H T+T- 


110-150 


4.7 |320j 




110-145 


4.9 [321] 


5.1 [319] 


H -> bb 


110-130 


4.7 |322j 




110-135 


5.0 |323] 


5.1 [319] 


H -> ZZ -> £+£-£'+£'- 


110-600 


4.8 |324j 


5.8 |3l7j 


110-600/160 


5.1 |325] 


5.3 |319] 


H -> W+W~ -> l + v£~v 


110-600 


4.7 |326j 


5.8 |3l7j 


110-600/160 


4.9 |327] 


5.1 [319] 


H ->■ ZZ -> l+l-TT 








190-600 


4.7 [328] 




H -»■ ZZ 


200-600 


4.7 |329j 




130-163,200-600 


4.6 |330] 




H ZZ -»■ £+£-z/z/ 


200-600 


4.7 |331j 




250-600 


4.6 [332] 




H W+W- ^i/gg 


300-600 


4.7 |333j 











150 GeV < M H < 200 GeV, and (c) 200 GeV < M H < 600 GeV. 

The mass range above 600 GeV, which is highly disfavoured by electroweak precision measurements, 
was not considered, because the theoretical tools or recipes described in Refs. [25l[26l|5ll[55] (see also 
references therein) for a proper description of the extremely broad Higgs-boson lineshape and their 
interference with irreducible background for such large masses were not available to the experiments 
until June 2012. 

Both experiments published their findings based on the complete 2011 data set |314[I315] . These 
data were partially reanalyzed with better reconstruction and identification tools, and optimized analysis 
techniques, resulting in a significantly increased sensitivity. In combination with the analysis of the data 
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set collected in 2012 up to June, the reanalyzed 2011 data are published in Refs. |317[I319] . CMS only 
considered the mass region below 160 GeV in the 2012 data analysis and reanalyzed the 2011 data in 
the channels H — > 77, H — > bb, H — > t + t~ , and H — > LL — > £ + £~£' + £'~. ATLAS reanalyzed the channels 
H — > 77 and H — > LL — > £ + £~£' + £'~ in the 2011 data, and so far has only investigated those two channels 
and the H — > W + W~ —> e/iz/z/ decay mode in 2012 data. 

In the low-mass region 110— 150 GeV the decays to photons, to a pair of tau leptons or bottom 
quarks, and to a pair of weak gauge bosons (WW and ZZ), which decay to four leptons, were considered 
by both experiments. The decays H — > 77 and H — > ZZ — > £ + £~£' + £'~ , with very low rate, provide an 
excellent resolution for the reconstruction of the invariant mass of the decay products of the Higgs- 
boson candidate with an accuracy of about 1—2%. In the decay mode H — > W + W~ — > £ + v£~v, which 
has a comparably high production rate, only the transverse mass can be reconstructed. The decays 
H — > t + t~ and H — > bb allow for reconstructing the mass of the Higgs-boson candidate with limited 
resolution in the range of 10—30%, due to the presence of several neutrinos in the final state or the 
limited energy resolution for the reconstruction of jets compared to charged leptons and photons. CMS 
also considered the decay H — > ZZ — > l + £~qq for mass values larger than 130 GeV, where the mass can 
be reconstructed with an accuracy of roughly 10%. 

In the intermediate-mass region 150—200 GeV the sensitivity of the search is dominated by the decay 
mode H — > W + W~ — > £ + v£~v, which has a large production rate, but allows only for reconstructing the 
transverse mass. The search for H — > LL — > £ + £~£' + £'~ with low rate, but excellent mass reconstruction 
capabilities provides a sensitivity that is roughly one order of magnitude smaller. 

In the high- mass region 200—600 GeV various final states arising from the decay of the Higgs boson 
into a pair of electroweak gauge bosons are considered. In addition to the final states that domi- 
nate the sensitivity in the intermediate- mass region, the decays H — > ZZ — > £ + £~qq, H — )■ ZZ — >■ £ + £~vv 
are searched for by both experiments as well, while H — > LL — > £ + £~tt is analyzed only by 
CMS, and H W+W" ->• £vqq only by ATLAS. The signal rates for H ->■ ZZ ->• £+£'£' +f~ and 
H —7- W + W~ — > £ + v£~v get quite small for large Mh, and hence the sensitivity is increased by in- 
cluding decay modes of the weak gauge bosons with larger branching ratios. These final states are not 
considered for smaller masses as the backgrounds, dominated by W and Z production in association 
with jets, increase in importance in this region and cannot be suppressed sufficiently or cannot be 
estimated with sufficient precision from control regions in data to the required level. Additionally the 
width of the reconstructed Higgs-boson mass distribution is no longer dominated by the experimental 
resolution, but by the natural width in the H — > ZZ — > £ + £r£' + £'~ channels for masses above 350 GeV. 
Hence the relative sensitivity to the other decay modes, which have a worse mass resolution or only 
allow an approximate mass observable to be reconstructed, gets weaker. 

The channels H — > W + W~ — > £ + v£^v and H — > LL — > £ + £~£' + £'~ contribute significantly to the over- 
all sensitivity over the full mass range. For small values of Mh it is important to lower the threshold 
for the transverse momenta of the selected leptons as much as possible, limited often by trigger rate 
requirements, and to retain simultaneously a low misidentification rate. 

• In the H — > W + W~ — > £ + v£~v analysis an excess of events with two leptons of opposite charge and 
large missing transverse momentum is searched for. Events are divided into separate categories 
according to jet multiplicity and lepton flavours, where the two-jet category has selection criteria 
designed to enhance sensitivity to the VBF production process. ATLAS so far has only considered 
the e/i final state for the 2012 data analysis due to the higher pile- up conditions with worse 
resolution for the missing transverse momentum and the correspondingly increased background 
from Drell-Yan production in the same-lepton-flavour final state. In events with no jets, the 
main background stems from non-resonant WW production; in events with one jet, the dominant 
backgrounds are from WW and top-quark production, and in events with two jets from top-quark 
production. In CMS in the analysis of the 2011 data MVA classifiers are trained for a number of 
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Figure 24: Distributions of the final discriminant in the H — > W + W — > £ + u£ v search at 8 TeV in 2012 
in ATLAS [SH] (left) and CMS [HH] (right). 



Higgs-boson masses to improve the separation of signal from backgrounds, and a search is made 
for an excess of events in the output distributions of the classifiers. In the 2012 analysis CMS uses 
sequential cuts, and the di-lepton invariant mass is used as final discriminant. ATLAS performs a 
cut-based analysis, optimized for different Mr ranges, and the WW transverse-mass distribution 
is used as the final discriminant. All background rates, except for very small contributions from 
WZ, ZZ, and W7* are evaluated from data. The distribution of the final discriminants in the 
analysis of the 8 TeV data show a clear broad excess in ATLAS and a less significant broad one 
in CMS (see Figure M>- 



• In the H — > ZZ — > £ + £~£' + £'~ channel an excess in the invariant-mass spectrum of the four iso- 
lated leptons over a small continuum background is searched for. The 4e, 4/i, 2e2/i final states 
are analyzed separately, since there are differences in the four-lepton mass resolutions and the 
background rates arising from jets misidentified as leptons. The main irreducible background 
from non-resonant ZZ production is mainly estimated from simulation. The smaller reducible 
backgrounds from Z+jets production, which mostly impacts the low four-lepton invariant-mass 
region, and top-quark pair production are estimated from control regions in data, or at least 
the normalization is validated in dedicated control regions. The final discriminant in ATLAS 
is the four-lepton invariant-mass spectrum. CMS exploits the different kinematics and angular 
correlations to discriminate further the signal process from the irreducible ZZ background. The 
information from the five angles, which describe the production and decay kinematics in the 
Higgs-boson candidate rest frame, are combined in the so-called MELA observable |334| . CMS 
uses a two-dimensional final discriminant consisting of the four-lepton invariant mass and the 
MELA output {Kd). A clear excess in the four-lepton invariant-mass spectra is observed in both 
experiments at masses of 125 GeV in ATLAS and 126 GeV in CMS, respectively (see Figure 125]) . 



In addition to the two channels discussed above, where the full mass range is covered, the channel 
H — > ZZ — > £ + £~vv provides the largest sensitivity in the high-mass range and dominates the sensitivity 
for Mh values larger than 300 GeV. 
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Figure 25: Four-lepton invariant-mass spectra in ATLAS |317j (left) and CMS [319J (right) for the 
combination of the 2011 and 2012 data. 



• In the search for H — > ZZ — > £ + £~vv events with a pair of electrons or muons with opposite 
charge, consistent with the decay of an on-shell Z, and large missing transverse momentum are 
selected. After applying additional kinematical requirements, which depend on the Higgs-boson 
mass hypothesis, a transverse-mass observable is used as the final discriminant. The background 
predictions for weak-gauge-boson pair production are obtained from simulation, whereas the other 
backgrounds are determined from control samples in collision data. A significant fraction of the 
selected signal events stems from the decay H — > W + W~ — > £ + u£~u especially for lower masses. 
In order to avoid double counting of signal yields, the selection is designed to have no overlap with 
the dedicated H — > W + W~ — > £ + v£rv analysis. The normalization of most background processes 
is estimated from control regions in data. 

• In the search for H — > ZZ — > £ + £~qq events with a pair of electron or muons and a pair of jets, 
which each have an invariant mass consistent with Mz, are selected. The events are categorized 
according to the number of b-tagged jets in the final state (none or at least one in ATLAS; 0, 1, or 
2 in CMS), in order to profit from the relatively large rate of b-jets from Z-boson decays present in 
the signal compared to the rate of b-jets found in the Z+jets background. Additional kinematical 
requirements are applied depending on the Higgs-boson mass hypothesis. CMS also exploits the 
different angular correlations between the final-state objects for the signal process compared to the 
background processes. The final discriminating variable is given by the invariant mass of the two 
leptons and two jets. CMS also searches in the mass range below 2Mz, where the mass- window 
criteria on the di-lepton mass is replaced by an upper bound. The dominant background from 
Z+jet production is estimated from control samples in data, defined by sidebands in the di-jet 
invariant mass. 

• The search for H — > W + W~ — > £vqq (only ATLAS) requires events with an isolated lepton and 
missing transverse momentum and two jets with an invariant mass compatible with a W boson. 
The mass of the Higgs-boson candidate can be constructed by solving a quadratic equation to 
solve for the component of the neutrino momentum along the beam axis exploiting the W mass 
constraint for the lepton plus neutrino system. The analysis further classifies events by lepton 
flavour and by the number of additional jets (0, 1, or 2), where the two-jet channel is optimized 
for the VBF production process. The final discriminant is given by invariant mass of the £^qq 
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Figure 26: Distributions of the final discriminants in the ATLAS searches [335] in the high-Mr^ re- 
gion: H -»■ ZZ -»■ £ + £"z/z/ (top-left), H -»■ ZZ -> (top-right), H ZZ -»■ £+£"bb (bottom-left), 
and H — )■ W + W~ — > tvqq (bottom-right). 



system. The background is modeled with a smooth function. 

• The search for H — > ZZ — > £ + £~tt (only CMS) selects events with electron or muon pairs, which 
are consistent with stemming from the decay of an on-shell Z boson, and a pair of tau-lepton candi- 
dates, which either decay hadronically or into an electron or muon. Final states with only electrons 
or muons are not considered, as they are already included in the dedicated H — > ZZ — > £ + £~£' + £'~ 
search. The final discriminant is given by the distribution of the di-lepton-di-tau invariant mass, 
constructed from the visible products of the tau-lepton decays, neglecting the effect of the accom- 
panying neutrinos. The mass resolution for Higgs-boson candidates is 10—15%. The dominant 
background from ZZ production is estimated from simulation. Other backgrounds are estimated 
from control regions in data. 

The final mass distributions (examples from ATLAS are shown in Figure 126]) agree with the expectations 
from background processes. 

In the low- mass range in addition to H — y W + W~ — y £ + v£~v and H — y ZZ — y £ + £~~£ l+ £'~ the search 
is performed for H — y 77, H — y t + t~ , and H — y bb. 

• The H — y 77 analysis selects events with two isolated photons with large transverse momenta. The 
dominant reducible backgrounds are multi-jet production and single-photon production in associ- 
ation with jets with cross sections several orders of magnitude larger than the signal cross section. 
These backgrounds can be suppressed to a level of less than 25% due to the excellent discrimi- 
nation capabilities of the LHC detectors between photon and jets. The irreducible background 
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Figure 27: Di-photon invariant-mass spectra after the inclusive selection in ATLAS |317] (left) and 
CMS [SSI (right) for the combination of the 2011 and 2012 data. 



from di-photon production can only be discriminated from the signal process by an excellent re- 
construction of the invariant mass of the di-photon system, which depends on the measurement 
of the energies of the photons and their opening angle. In both experiments the contribution 
from the resolution in the measurement of the opening angle was found to be negligible. The 
di-photon mass spectra after requiring two isolated photon candidates for the data from 2011 
and 2012 are shown in Figure [271 The sensitivity of the search can be enhanced by splitting the 
event sample in various categories with different ratios of expected signal over background event 
yields and different signal mass resolutions. Both experiments consider event classes exploiting 
the typical VBF signature of two additional jets with large separation in pseudorapidity and large 
di-jet invariant mass. ATLAS uses nine additional categories based on the pseudorapidity of each 
photon, whether it is reconstructed as a converted or unconverted photon, and the momentum 
component of the di-photon system transverse to the di-photon thrust axis. CMS uses four addi- 
tional categories depending on the score of a boosted decision tree. The background in the signal 
region is estimated from a fit to the observed di-photon mass distribution in data. A clear excess 
in the invariant di-photon mass spectra is observed in both experiments at masses of 126 GeV in 
ATLAS and 125 GeV in CMS, respectively. 

• The H — > t + t~ analysis searches for a broad excess in the reconstructed di-tau invariant- mass 
distribution, where the mass reconstruction is mostly based on the ideas in Refs. [33611337] with a 
resolution between 10—30%. Events are classified according to the tau-lepton decay modes (into 
an electron, a muon, or hadrons plus neutrinos). ATLAS considers all decay- mode combinations, 
whereas CMS has not investigated yet the double electron and double hadronic final state. De- 
pending on the final state the event samples are further divided into exclusive subcategories, which 
are optimized for the sensitivity to a particular production mode. The category optimized for the 
VBF signature provides the largest sensitivity. Events failing these selection requirements are 
separated further, e.g., in a category with two jets optimized for associated production HW(HZ), 
with one jet with large transverse momentum for associated production and production in gluon 
fusion, and with either no jets or with one with a small transverse momentum aiming at the 
gluon-fusion production process. The main irreducible background arises from Z — > tt produc- 
tion, whose di-tau invariant-mass distribution is derived from data by selecting Z — > //// events, 
in which the reconstructed muons are replaced with reconstructed particles from the decay of 
simulated r leptons of the same momenta. The reducible backgrounds from W+jets, Drell-Yan, 
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Figure 28: Distributions of the final discriminant for the search in H — > t + t combined for 2011 and 
2012 data (left) and H -» bb for 2012 data (right) in CMS [3T9] . 



and multi-jet production are also evaluated from control samples in data. The distribution of 
the di-tau mass distribution for 2011 and 2012 data in CMS is shown in Figure |2"51 (left). No 
significant deviation from the background expectation is observed. 

• The H — > bb analysis selects events produced in association with a W or Z decaying via Z — > £ + £~ , 
Z — > vv, or W — > tv. The selection requires two b-tagged jets and either two leptons, one lepton 
and missing transverse momentum, or large missing transverse momentum. CMS selects events 
with a high transverse momentum of the di-jet system in order to improve the mass resolution 
and to suppress background processes. A multivariate classifier is trained for different Mh values 
and its output is used as the final discriminant in CMS (see Figure |28| right). ATLAS subdivides 
the events in categories according to the transverse momentum of the weak gauge boson and uses 
the di-jet invariant mass of the b-tagged jets as the final discriminant. No significant deviation 
from the background expectation is observed. 

8.3 Excluded mass ranges 

In most search channels and over almost the full investigated mass range from 110—600 GeV no signif- 
icant excess with respect to the background-only hypothesis is observed. This allows for the exclusion 
of a large part in the mass range of the SM Higgs boson. In the combination of the different search 
channels the ratio of the cross section for the different production modes and the ratio of the decay 
branching ratios are taken from the SM prediction. Upper limits on the signal strength /i, which relates 
the excluded cross section to that for Higgs-boson production in the SM, are derived based on the pro- 
file likelihood ratio |109] as test statistic using the CLs technique |106] to avoid exclusion of parameter 
space where no sensitivity is expected. Correlations in the theoretical and systematic uncertainties, 
which are included via nuisance parameters in the likelihood function, among the various search chan- 
nels and event categories are taken into account (see Refs. |114|.l314p315| for details). The exclusion 
limits obtained in both experiments in the various search channels and their combination from the 
first analysis of the full 2011 data set are shown in Figure |29| for the combination in each experiment 
from the (re-)analysis of 2011 and 2012 data in Figure |30| and both are summarized in Table ITOl 
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m H l GeV l Higgs boson mass (GeV) 



Figure 29: Exclusion limits on the signal-strength parameter /x as a function of the hypothetical Higgs- 
boson mass in the individual channels and their combination obtained from the analysis of 2011 data 
by ATLAS 032] (left) and CMS [215] (right). 




Figure 30: Exclusion limits on the signal-strength parameter // as a function of the hypothetical Higgs- 
boson mass obtained from the analysis of 2011 and 2012 data by ATLAS [317] (left) and confidence 
level CLg for rejecting the signal+background hypothesis as a function of the hypothetical Higgs-boson 
mass by CMS (right). 



Already with the data set from 2011 a large mass range is excluded. However, in the low- mass range 
the observed limit is weaker than the expected one, which extends from 118— 560 GeV (by applying a 
simple "or" of the two experiments). After adding the 2012 data set and re-analyzing the 2011 data 
the expected exclusion sensitivity extends over the mass range of 110—582 GeV. The mass range that 
is excluded neither by ATLAS nor by CMS, is the interval from 122— 127 GeV. The large deviations of 
the observed from the expected limits are caused by the observed excesses of events discussed above. 
For a large mass range, reduced production cross sections of 20% of the SM prediction are excluded as 
well, and the SM Higgs-boson hypothesis can be excluded with very high confidence level. 
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Table 10: Mass ranges for the SM Higgs boson excluded by the combination of all searches of ATLAS and 
CMS: observed excluded range, expected excluded range for 2011 data alone {314(1315] , and combination 
of reanalyzed 2011 data and 2012 data [51711515] . 





2011 df 
excl. M H [GeV] 


it a 

exp. excl. M H [GeV] 


2011+2 
excl. Mn [GeV] 


012 data 

exp. excl. M H [GeV] 


ATLAS 


111.4-116.6, 119.4-122.1 
129.2-541 


120-560 


111-122, 131-559 


110-582 


CMS 


127-600 


118-543 


110-121.5, 127-600 


110-543 




Figure 31: Local p-values as function of the hypothetical Higgs-boson mass in ATLAS [317] (left) and 
CMS [319] (right). 

8.4 Observation of a new particle 

When testing the compatibility with the background-only hypothesis, already in the 2011 data, tanta- 
lizing hints for the possible production of a new particle were observed, which corresponded to local 
significance of 2.9 at a mass of 126 GeV in ATLAS [335] and a local significance of 3.1 at a mass of 
124 GeV in CMS |315j . These local significances are already larger than the ones shown by the individual 
Tevatron experiments and quite similar to the combined significance published by CDF and DO. When 
taking into account the "look-elsewhere effect" in the whole mass range considered from 110— 600 GeV 
the global significances were found to be roughly 1 in ATLAS and 1.5 in CMS. After adding the 2012 
data available up to June to the analysis the local significances (p-values) increased (decreased) to 5.9 
(1.7 x lO" 9 ) at 126.5 GeV in ATLAS and 5.0 (2.8 x 10" 7 ) in CMS (see Figure EH). An overview of the 
observed and expected significances in the various search channels and their combination is given in 
Table [Til ordered by the size of the excess. The observed significance is due to the excesses observed 
in the channels with high mass resolution H — > 77, H — > ZZ — > £ + £~£' + £'~ and with large production 
rate H — > W + W~ — > £ + v£~v. No hints for the decay into a pair of fermions are observed yet, though 
the expected significance in the combination of H — > t + t~ and H — > bb is equal to the expected one 
in H — > W + W~ i^v^v in the CMS analysis. Assuming that the observed signal is due to the SM 
Higgs boson, ATLAS observes an upwards fluctuation of roughly one standard deviation compared to 
the expectation at a mass of 126.5 GeV and CMS observes a downwards fluctuation of roughly one 
standard deviation at a mass of 125.5 GeV. The global significance is 5.1 (5.3) in ATLAS for the mass 
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Table 11: Local significances in different search channels in ATLAS [317J and CMS |319j . Z denotes the 
maximum significance observed at the mass value M max z when scanning the hypothetical Higgs-boson 
mass in the hypothesis test, and Z exp . is the median expected significance for a SM Higgs boson at the 
same mass value. An entry "-" means that the information is not publicly available. 
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Figure 32: Best fit values for the signal strength jj as function of the hypothetical Higgs-boson mass in 
ATLAS §H\ (left) and CMS [3l9] (right). 

range 110-600 (110-150) GeV and 4.5 (4.6) in CMS for the mass range 110-145 (110-130) GeV. These 
findings encouraged the two experiments ATLAS and CMS to proclaim the discovery of a new particle 
and to entitle their publications [317jl319j with "Observation of a new particle / boson" . 

8.5 Anatomy of the new particle 

From the observation of the decay into a pair of particles with identical spin and with vanishing sum 
of the electric charges it can be concluded that the particle is electrically neutral, has integer spin, and 
hence is a neutral boson. The observation of the decay in two massless identical vector bosons, H — > 77, 
excludes the hypothesis of a spin-1 particle according to the Landau- Yang theorem [338J. 

A first crucial test of the compatibility of the observed event yields with the prediction for the SM 
Higgs boson is provided by a fit of the signal-strength parameter fi. The overall signal strength relative 
to the SM prediction fi, which assumes the ratio of production cross sections and ratio of branching 
ratios as predicted by the SM, is shown in Figure 1321 as function of the hypothetical Higgs-boson mass. 
In both experiments the largest signal strengths of 1.4 ± 0.3 observed at a mass of 126.0 GeV (ATLAS) 
and of 0.87 ± 0.23 observed at a mass of 125.5 GeV (CMS) are compatible with unity. In a next step 
the best signal strength is determined for each decay mode of the new boson, which only assumes 
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Figure 33: Best fit values for the signal strength \i in the different decay modes as function of the 
hypothetical Higgs-boson mass in ATLAS j3T7] (left) and CMS jSTTO] (right). 

that the different production modes contribute as predicted by the SM. The fit results for the different 
decay modes, where the hypothetical Higgs-boson mass is fixed to the value that yields the largest 
overall signal strength, are also consistent with the SM prediction (Figure |33|) within the still large 
uncertainties. 

The compatibility of a simultaneous determination of the best signal strength /i and the best mass 
value Mh in several final states is demonstrated in Figure EH The best mass value Mh is determined 
in both experiments from the observed mass spectra in the H — > 77 and H — > ZZ — > £ + £~£' + £'~ final 
states. The signal strengths in each final state (decay modes and classification categories) are treated 
as independent nuisance parameters. The results for the mass of the new boson, as obtained by the two 
collaborations are 

ATLAS: 126.0 ± 0.4(stat.) ± 0.4(sys.) GeV [317] and CMS: 125.3 ± 0.4(stat.) ± 0.5(sys.) GeV [3T9] . 

In both experiments the systematic uncertainty is dominated by the knowledge of the absolute energy 
scale of photons and to a lesser extent of that of electrons. The mass values determined in the different 
final states within one experiment as well as the mass values determined in the two experiments are 
consistent with each other. Already now the precision of the determination of the mass is better than 
one percent. Note that the measured mass value is perfectly compatible with the Mh range preferred by 
the fit of the SM to electroweak precision measurements (see Section 12.51) . The particle discovery, thus, 
is not only a triumph of the interplay of LHC measurements and corresponding theoretical predictions, 
but also a spectacular success of electroweak precision physics — assuming the new particle is in fact the 
Higgs boson of the SM. 

9 Conclusions and outlook 

Almost 50 years after the postulate of an elementary scalar boson, which in the standard theory of 
electroweak interactions — the Glashow-Salam- Weinberg model — is the companion of the mechanism 
to describe the spontaneous breaking of the electroweak symmetry, a new neutral boson with a mass of 
126 GeV has been discovered in the search for the SM Higgs boson at the LHC. So far all observed event 
yields in the different final states are consistent with the expectations for the SM Higgs boson, but also 
with predictions for a Higgs boson in several extended models. Assuming that the excess observed at 
Tevatron in the bb final state is not due to a statistical fluctuation, but due to the production of the 
particle discovered at the LHC, the decay of the new particle into a pair of b-quarks and probably also 
into a pair of tau leptons should be observed at LHC in the near future. 
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Figure 34: Best fit values for the signal strength /i versus hypothetical Higgs-boson mass Me in 
ATLAS [3H] (left) and CMS [319] (right). 



The allowed mass range for the Higgs boson could be confined unambiguously by the tremen- 
dous and challenging efforts in the searches at the colliders LEP, Tevatron, and LHC during the last 
decades: The combination of ALEPH, DELPHI, L3, and OPAL results exclude M H = 0-114.3 GeV, 
the preliminary combination of CDF and DO results additionally exclude Mh = 147— 180 GeV (and 
Mh = 100—103 GeV), and finally at least one of ATLAS and CMS excludes the SM Higgs boson in the 
ranges Mh = 110— 122 GeV and Mh = 127— 600 GeV. In significant parts of the excluded mass ranges 
a Higgs boson with reduced couplings with respect to the SM prediction is also excluded. 

The role of the Higgs boson in the zoo of elementary particles is rather special. It is the only 
elementary spin-0 boson and shows a very distinctive interaction pattern with all other particles. The 
scalar nature of the new boson can be affirmed by excluding the alternative integer-spin hypotheses 
from the measurement of angular correlations in production and decay of the new boson. Assuming it 
is a spin-0 boson this will probably be possible with the data collected at LHC until the end of 2012. 
In order to verify the consistency with being a Higgs boson, the Lorentz structure of the interactions 
with all particles to which a coupling has been or will be observed, in particular to the massive gauge 
bosons, has to be determined. This includes the determination of its CP nature, which in the SM is 
predicted to be CP even. In models with extended Higgs sectors CP-odd states or states of mixed CP 
nature are often predicted. 

The idea of the Higgs mechanism interprets particle masses as an interaction strength of each 
elementary massive particle to the vacuum expectation value of the Higgs field, whose particle excitation 
represents the Higgs boson. As a result, all massive elementary particles couple to the Higgs boson 
proportional to their masses. In a first approach, the coupling strengths can be determined from the 
measurement of the products of production cross section times branching ratio, which can be extracted 
from the observed event yields in the different final states. Later, more sophisticated coupling analyses 
based on effective field theories, quantifying non-standard effects in a model-independent way, may 
be carried out with higher precision. For a mass of 126 GeV, for which the SM predicts the total 
Higgs-boson width to be 4MeV, the total width cannot be determined at the LHC with sufficient 
precision. As a consequence only ratios of couplings can be determined at the LHC in a completely 
model-independent way. However, with mild assumptions, e.g. that the couplings to massive weak 
gauge bosons are not larger than predicted by the SM (valid in many extended Higgs sectors) , absolute 
coupling strengths can be measured as well. At a future e + e~ collider absolute coupling values can be 
determined in a completely model-independent way and probably with a significantly better precision. 
Finally, a determination of the trilinear Higgs-boson self-coupling would be highly desirable in order to 
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at least partially reconstruct the Higgs potential. Owing to the low signal rates that are expected for 
Higgs-boson pair production this will be very challenging at the LHC, both in its run at a CM energy of 
14 TeV and in its high-luminosity phase, but also future e + e~ colliders will have to be pushed to their 
limits to this end. 

Assuming the observed particle will be identified as a Higgs boson, the power to discriminate the 
SM Higgs boson from Higgs bosons in extended models depends on the achievable precision in the 
coupling determination. The possibly observed qualitative or quantitative differences between SM 
theory predictions and experiment in Higgs-boson observables could indicate the way, how the SM may 
have to be extended, and tell which models are clearly disfavoured or can be ruled out. In addition, 
the future research programme includes the search for additional Higgs bosons in non-standard scalar 
sectors, the search for other new particles, and the investigation of processes like weak-gauge-boson 
scattering that are most sensitive to the mechanism of electroweak symmetry breaking. 

Any model for the unification of the fundamental strong and electroweak forces will predict new 
heavy particles, which in turn receive their mass most likely via spontaneous symmetry breaking. 
Accepting the Higgs mechanism as the correct modelling of electroweak symmetry breaking, it is, thus, 
reasonable to expect that new Higgs bosons arise in SM extensions. Conversely, precision measurements 
of the or more Higgs bosons serve as a window to non-standard theories or even grand unification. 
Supersymmetric theories are typical examples sharing this feature. 

Apart from these general considerations, there are at least two known phenomena where our estab- 
lished theoretical world, which consists of the Standard Model of particle physics and general relativity, 
finds its limitations: The unknown constitution of Dark Matter and the enigmatical role and size of 
Dark Energy in the universe. Clarifying the nature and profile of the new particle is of utmost im- 
portance especially in the light of those conundrums, since extended Higgs models can accommodate 
Dark Matter candidates and Higgs fields could in fact play a key role in the inflationary phase of the 
cosmological expansion that is driven by Dark Energy. 

Particle physics has reached a crossroad, and Higgs precision physics will probably contribute a lot 
to steer the field in future directions. 
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